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CHAPTER  I 

BACKGROUND  AND  MOTIVATION 


The  Defense  Nuclear  Agency  and  the  importance  of  arms  control 
are  discussed.  Passive  radiation  detection,  including  neutron  and 
gamma  signatures,  is  described.  The  report  is  outlined. 


1 . 1  MOTIVATION 

One  of  the  Defense  Nuclear  Agency's  (DNA)  missions  is  to 
develop  and  evaluate  new  technologies  that  could  be  applied  to  the 
arms  control  verification  process.^  The  proliferation  of  nuclear 
weapons  technology,  especially  to  the  Third  World,  and  the  ongoing 
reduction  in  the  nuclear  stockpile  of  the  former  Soviet  Union  as 
part  of  the  Strategic  Arms  Reduction  Treaty  (START)  have  made  arms 
control  an  extremely  important  issue.  Now,  more  than  ever  it  would 
be  useful  to  develop  reliable  and  non- intrusive  technology  that 
could  be  used  in  the  field  by  inspection  teams  to  evaluate  whether 
the  number  of  nuclear  weapons  are  within  the  specific  guidelines 
delineated  by  the  nuclear  forces  reduction  treaties,  specifically 
START.  Inspectors  must  be  able  to  determine  both  the  presence  of 
strategic  offensive  arms  and  whether  those  arms  have  no  more  than 
the  maximum  allowable  number  of  warheads.  Other  potential 
inspection  technologies  include  radiographic,  gravimetric,  and 
acoustic  processes,  among  others^.  Although  these  methods  may  be 
effective,  they  can  be  relatively  expensive  and  may  often  reveal 
sensitive  weapon  design  information.  This  is  a  problem  because 
countries  do  not  want  to  divulge  unique  weapon  design 
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specifications  when  not  obligated  to  do  so  by  the  treaty. 
Furthermore,  the  concept  of  reciprocity  does  apply  to  the 
verification  techniques  allowed  under  the  current  nuclear  forces 
reduction  treaties. 

1.2  PASSIVE  RADIATION  DETECTORS 

DNA  currently  believes  that  passive  radiation  detectors  might 
non- intrusively  provide  necessary  information  for  nuclear  treaty 
verifications.^  In  July  1989,  this  concept  was  evaluated  outside 
the  U.S.  Government  in  a  series  of  experiments  conducted  under  the 
auspices  of  the  Natural  Resources  Defense  Council  (NRDC)  and  the 
Academy  of  Sciences  of  the  U.S.S.R. The  Soviet  Navy  provided  the 
cruiser  SlsLVSt  armed  with  a  single  nuclear —armed  SS— N— 12  sea- 
launched  cruise  missile  (SLCM) .  Using  numerous  gamma-ray 
detectors,  the  investigators  detected  the  presence  of  a  nuclear 
warhead  without  the  need  to  have  physical  access  to  the  weapon. 

The  Naval  Research  Laboratory  (NRL)  has  proposed  the  use  of 
neutron  signatures  induced  as  a  result  of  the  bombardment  of  the 
earth  by  cosmic  rays  for  Nuclear  Treaty  verification  applications.^ 
The  signature  of  the  neutrons  emitted  by  a  body  is  potentially 
unique  and  depends  on  the  internal  structure  and  materials  of  the 
object.  As  a  result,  induced  neutron  measuring  devices  might  be 
applied  to  verifications  processes.  Currently  NRL  is  studying  the 
use  of  a  high  efficiency  Helium-3  neutron  detector  to  measure  these 
signatures  accurately.® 

DNA  now  believes  direct  measurements  of  weapon  neutron 
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signatures  offer  more  promise  for  verification  than  the 
aforementioned  induced  radiation  detection  techniques.^  Nuclear 
warheads,  because  of  their  composition  of  enriched  uranium  and 
plutonium,  produce  a  characteristic  neutron  signature.  An  accurate 
measurement  of  this  signature  should  yield  information  concerning 
the  presence  of  nuclear  warheads. 

1.3  POSSIBLE  NEUTRON  DETECTORS 

Neutron  detection  is  inherently  difficult  and  complex  owing  to 
the  physical  characteristics  of  neutron  radiation.  Neutrons  are 
neutral  particles  that  interact  through  a  variety  of  mechanisms  to 
form  secondary  charged  particles.  Neutron  detection  relies  upon 
the  measurement  of  these  secondary  particles  in  order  to  deduce  the 
amount  and  energy  of  the  original  neutron  radiation.*  However,  the 
probability  of  a  given  nuclear  interaction  is  dependent  upon  the 
kinetic  energy  of  the  incident  neutron.  Neutron  capture  dominates 
at  lower  energies  and  causes  alpha,  proton,  or  gamma  emissions, 
while  elastic  scattering,  which  generates  heavy  recoil  ions,  is 
predominant  among  higher  energy  neutrons.  Because  typical  neutron 
fields  of  interest  cover  a  broad  spectrum  of  energy,  from  0.025  eV 
(thermal)  to  about  6  MeV,  accurate  neutron  detection  can  be  quite 
complicated  due  to  the  mixture  of  mechanisms  involved. 

There  are  several  types  of  commercially  available  neutron 
detectors,  each  having  its  own  unique  advantages  and  drawbacks.  Two 
of  the  devices  that  have  been  identified  as  possible  candidates  for 
development  as  field  instruments  are  the  Helium-3  (He-3) 
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proportional  counter  and  the  neutron  bubble  dosimeter.  The  He-3 
detector  indirectly  measures  the  presence  of  neutrons  using  the 
ionizing  gas  Helium-3.  When  a  neutron  reacts  with  He-3  gas,  the 
following  reaction  occurs: 

jHe^*  +  on^  -  iH"  +  ip^  (1-1) 

The  products  of  the  reaction  are  a  triton  (i.e.,  tritium  atom  with 
an  extra  electron) and  a  proton.  When  a  potential  is  applied  across 
a  He-3  filled  chamber,  each  neutron  interaction  produces  a  pulse. 
He-3  gas  has  an  extremely  high  thermal  neutron  capture  cross- 
section,  but  a  small  fast  neutron  capture  cross  section.  As  a 
result,  the  device  must  be  equipped  with  a  moderator  when  used  to 
detect  a  multi-energy  neutron  spectrum  so  that  the  response  of  the 
detector  will  not  be  severely  limited. 

The  basic  operation  of  neutron  bubble  detectors  is 
fundamentally  the  same  as  that  of  bubble  chambers,  the  only 
difference  being  that  instead  of  filling  a  chamber  completely  with 
superheated  liquid,  bubble  detectors  contain  many  superheated 
liquid  droplets  which  act  as  individual  microscopic  bubble 
chambers.  Both  devices  operate  on  the  principle  first  explained  by 
Glaser:’  if  enough  energy  is  transferred  from  a  moving  ion  to  a 
superheated  liquid,  then  complete  vaporization  of  the  liquid  is 
possible.  The  recently  developed  bubble  detectors  utilize 
thousands  of  independent  superheated  liquid  droplets  suspended  in 
a  semi-viscous  matrix.  When  neutrons  interact  with  the  droplet 
material,  the  secondary  ions  deposit  sufficient  energy  to  nucleate 
them  into  macroscopically  observable  vapor  bubbles.  The  number  of 
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bubbles  can  be  directly  related  to  the  amount  of  neutron  exposure, 
thus  providing  a  portable,  self-reading,  isotropic  detector  which 
is  simple  and  fairly  inexpensive*  By  varying  the  superheated 
liquid  material  or  the  amount  of  droplets  suspended  in  the  gel 
matrix,  different  energy  thresholds  and  sensitivities  can  be 
obtained.  The  Navy  has  conducted  extensive  research  evaluating  the 
bubble  detector  for  personal  dosimetry  applications.* 

There  are  currently  two  commercially  available  bubble 
dosimeters,  each  having  a  different  scheme  for  counting  the 
vaporizations.  The  Apfel  detector,*  a  product  of  Apfel 
Enterprises,  acoustically  monitors  the  device  during  neutron 
irradiation  and  electronically  records  each  event  by  means  of  a 
piezoelectric  transducer. “  Ing's  device, manufactured  by  Bubble 
Technologies  Incorporated  (BTI) ,  is  much  simpler.  The  bubbles  form 
in  a  more  rigid  matrix  and  are  held  in  place  after  vaporization  so 
that  they  can  be  visually  counted  either  by  eye  or  with  a  computer 
driven  optical  reader.  Unlike  the  Apfel  device,  the  BTI  detector 
can  be  recompressed  after  reading  to  condense  the  bubbles,  making 
it  reusable  and  thus  cheaper.  Both  versions  of  bubble  dosimeters 
have  been  shown  to  be  insensitive  to  gammas  below  6  MeV.^^ 

Although  the  He-3  proportional  detector  and  the  bubble  neutron 
detector  are  both  capable  of  counting  neutrons,  the  question  has 
arisen  as  to  which  device  is  better  able  to  carry  out  the  mission 
of  nuclear  treaty  verification.  In  order  to  fulfill  the  needs  of 
DNA,  the  neutron  detector  must  be  rugged,  easily  transportable,  and 
non- intrusive.  Detector  weight  becomes  an  extremely  important 
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issue  when  comparing  the  weights  of  current  proportional  counters 
to  the  weight  of  the  bubble  detector.  The  proportional  detectors 
being  evaluated  have  total  weights  ranging  from  four  to  35  pounds,* 
while  the  bubble  detector  weighs  only  0.125  pounds. 

Only  the  bubble  dosimeter,  specifically  the  BTI  device,  meets 
all  of  these  criteria.  Because  of  the  low  level  of  technology 
associated  with  the  device,  the  detector  requires  no  electronics. 
As  a  consequence,  the  device  is  much  smaller  and  more  reliable  in 
the  field  than  the  He-3  detector  or  the  Apfel  device.  Although  the 
electronics  associated  with  the  Apfel  device  are  minimal,  the 
dsvice  has  problems  with  background  noise  while  acoustically 
monitoring  bubble  formations.  Furthermore,  the  neutron  response 
of  the  bubble  detector  is  more  related  to  the  energy  of  the  source 
neutrons  than  the  He— 3  detector.  In  Fig.  1—1,  it  can  be  seen  that 
the  typical  response  for  a  BTI  bubble  detector  has  a  plateau  from 
about  0.3  MeV  up  to  about  2  MeV.  As  mentioned  above,  the  He-3 
detector  is  a  thermal  neutron  detector.  If  it  is  used  to  measure 
a  higher  energy  neutron  field,  then  it  must  be  equipped  with  a 
moderator  to  slow  the  fast  neutrons  down  into  the  thermal  range. 

Although  the  BTI  device  has  the  advantage  of  being  more 
portable  for  the  field  environment  than  the  He-3  detector,  the  He- 
3  detector  does  have  some  specific  advantages.  The  two  devices 
differ  drastically  in  their  means  of  counting  neutrons.  The  He-3 
detector,  because  of  its  design,  is  capable  of  achieving  much 
higher  counting  efficiencies  (provides  a  high  neutron  count  rate) 
than  the  BTI  device.  As  a  result,  the  counting  time  associated 
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FIGURE  1-1.  A  typical  response  as  a  function  of  neutron  energy 

for  a  neutron  bubble  dosimeter. 

with  statistically  accurate  results  is  much  smaller  for  the  He-3 
detector  than  for  the  BTI.  Because  of  its  higher  counting 
efficiency,  the  He-3  detector  would  be  better  suited  than  the  BTI 
device  for  the  detection  of  concealed  nuclear  weapons  where  the 
neutron  field  is  much  smaller.  However,  the  BTI  device,  because  of 
its  small  weight  and  low  level  of  technology,  still  remains  a 
potential  detector  for  verification  purposes. 

1 . 4  REPORT  OUTLINE 

Although  the  neutron  bubble  detector  has  been  classified  as  a 
possible  candidate  for  treaty  verification  applications,  several 


drawbacks  have  been  identified  with  the  use  of  bubble  dosimeters. 
First,  most  of  the  commercially  available  devices  have  a  strong 
temperature  dependence,  which  limits  their  use  to  constant 
temperature  environments.  As  seen  in  Fig.  1“2,  the  relative 
response  of  these  devices  increases  non^linearly  with  increasing 
temperature.  Second,  not  all  the  bubbles  produced  are  immediately 
visible,  but  instead  require  a  waiting  time  of  several  minutes  or 
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FIGURE  1-2.  The  relative  response  as  a  function  of  temperature 

for  an  uncompensated  BD-IOOR  neutron  bubble 
dosimeter. 

hours  depending  on  the  initial  droplet  size.  Third,  the  devices 
have  a  limited  dynamic  range  of  about  a  thousand  bubbles  per 
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exposure.  Because  of  the  statistical  nature  of  the  bubble 
formation  process  and  the  relatively  small  number  of  bubbles  often 
produced,  statistical  uncertainty  when  counting  also  needs  to  be 
considered. 


This  report  is  a  feasibility  study  for  the  Defense  Nuclear 
Agency  as  part  of  the  Arms  Control  Research,  Development,  Test,  and 
Evaluation  Program  (RDT&E)  and  will  evaluate  the  potential  of  the 
neutron  bubble  detector  for  START  verification  applications.  The 
report  will  focus  upon  the  effectiveness  of  the  device  as  a  neutron 
detector  and  will  attempt  either  to  solve  or  estimate  the  extent  of 
the  problems  previously  outlined. 

A  background  of  the  nuclear  reactions  that  occur  in  the 
detectors  is  provided  in  Chapter  II.  Using  a  theoretical  model 
developed  by  Dr.  Mark  J.  Harper,®  the  neutron^heavy  ion-*bubble 
process  is  presented. 

Chapter  III  investigates  the  use  of  different  droplet 
materials  as  a  solution  to  the  device's  temperature  sensitivity 
problem.  The  Harper  model  will  be  used  extensively  to  establish 
and  evaluate  possible  candidate  materials  theoretically. 

Several  experiments  were  conducted  at  the  National  Institute 

» 

of  Standards  and  Technology  (NIST)  to  validate  the  theoretical 
model.  The  setups,  results,  and  analyses  of  the  data  are  reviewed  <» 

in  Chapter  IV. 

Additional  experimentation  was  conducted  at  USNA  to  evaluate 
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further  the  problems  associated  with  the  device  and  its  potential 
as  a  reliable  neutron  detector.  Chapter  V  presents  the  setups, 

results,  and  analyses  of  the  data. 

Chapter  VI  draws  conclusions  concerning  the  feasibility  of 
using  a  neutron  bubble  detector  as  a  field  detector  for  nuclear 
treaty  verification  applications.  Possible  ideas  for  further 
research  are  also  discussed. 
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CHAPTER  II 


THEORETICAL  MODEL 

The  significant  neutron  reactions  in  typical  superheated 
liquid  droplet  materials  are  examined.  The  radiation- induced 
nucleation  of  bubbles  in  superheated  liquid  droplets  is  discussed. 
A  model,  developed  by  M.J.  Harper,*  which  accounts  for  all  probable 
neutronic  interactions  resulting  in  the  formation  of  bubbles  in 
superheated  liquid  droplets  is  outlined. 


2.1  NEUTRONIC  INTERACTIONS 

Several  conditions  must  be  met  for  neutrons  to  cause  a  bubble 
to  nucleate  from  a  superheated  liquid  droplet.  First,  some  type  of 
reaction  must  occur  between  the  neutron  and  the  liquid  material. 
The  probability  of  this  interaction  is  highly  dependent  upon  the 
kinetic  energy  of  the  incident  neutron,  the  overall  inclination  of 
the  target  atom  for  specific  types  of  interactions  (nuclear  cross 
section),  and  the  atomic  density  of  the  target  material.  Second, 
some  type  of  ionizing  particle  or  other  radiation  capable  of 
dissipating  energy  across  a  very  small  distance  must  be  created  as 
a  result  of  the  nuclear  reaction.  Neutrons,  being  electrically 
neutral,  are  not  capable  of  directly  ionizing  matter.  The 
secondary  radiations,  which  result  from  the  primary  neutron 
interactions,  deposit  their  energy  through  electronic  and  nuclear 
collisions  in  the  medium.  Third,  sufficient  energy  must  be 
deposited  within  a  small  enough  distance  for  a  permanent  phase 
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transformation  from  liquid  to  vapor  to  occur. 

Neutron  interactions  producing  secondary  ionizing  radiation 
are  classified  as  either  scattering  or  absorption.  Scattering 
reactions  produce  heavy  recoil  ions  (the  target  atom  stripped  of 
its  electrons)  and  a  deflected,  lower  energy  neutron.  Such 
reactions  are  either  elastic,  where  momentum  and  kinetic  energy  are 
conserved,  or  inelastic,  where  the  target  nucleus  is  left  in  an 
excited  state,  to  decay  later  by  beta  or  gamma  emission. 
Absorption  reactions,  however,  do  not  always  re-emit  a  neutron. 
The  energy  of  the  incident  neutron  may  appear  either  as  a  photon 
instantaneously  released  or  as  an  energetic  particle  ejected  from 
the  nucleus.  Fission  may  also  result  if  the  material  is  fissile. 
Spallation,  which  is  possible  at  energies  above  about  20  MeV, 
ejects  several  nuclear  fragments  and  particles. 

When  a  polyenergetic  neutron  radiation  field  characteristic  of 
a  nuclear  warhead  (0.025  eV  to  6  MeV)  is  incident  upon  a  neutron 
bubble  detector,  the  only  nuclear  reactions  possible  are  elastic 
scattering,  inelastic  scattering,  and  charged  particle  productions 
such  as  (n,p) ,  (n,He3)  ,  and  (n,Q!) . 

Refrigerants  are  typically  the  best  choice  for  the  superheated 
liquids  used  in  bubble  detectors.^^^  A  combination  of  low  boiling 
points  at  atmospheric  pressure,  high  heats  of  vaporization,  and 
favorable  physical  properties  make  refrigerants,  specifically 
chlorof luorocarbons,  fluorocarbons,  and  hydrocarbons,  ideal 
materials  for  radiation  induced  nucleation.  Since  carbon, 
chlorine,  fluorine,  and  hydrogen  are  the  only  elements  found  in 
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these  compounds,  the  probability  and  significance  of  all  reactions 
can  readily  be  determined  for  each  isotope.  If  a  reaction  is 
physically  possible  and  its  cross  section  is  non-trivial  for  the 
energy  range  considered,  then  the  effects  of  the  reaction  must  be 
considered. 

2 . 2  BUBBLE  NUCLEATION 

Several  theories  have  been  offered  as  to  the  exact  mechanism 
by  which  energy  converts  from  the  kinetic  in  charged  particles  to 
thermal  during  radiation-induced  nucleation.  However,  the 
currently  accepted  theory  is  that  of  the  "thermal  spike. 
According  to  Seitz,  the  originator  of  the  concept,  the  intense 
energy  deposition  along  the  charged  particle's  track  provides 
enough  localized  heating  in  the  superheated  liquid  to  cause 
nucleation. 

There  are  two  steps  in  bubble  formation  by  the  thermal  spike 
processes.  First,  a  small  "seed"  bubble  is  formed  when  a  small 
amount  of  superheated  liquid  is  directly  vaporized  by  the  energy 
transferred  from  the  recoil  ion.  Second,  the  bubble  grows  to 
visible  proportions.  Once  formed,  this  seed  bubble  begins  to 
expand  as  it  receives  more  thermal  energy  due  to  the  propagation  of 
the  thermal  spike's  shock  wave  through  the  surrounding  superheated 
liquid.  However,  the  bubble  must  expand  beyond  a  critical  size  if 
it  is  to  grow  and  eventually  consume  the  entire  liquid  droplet.  A 
vapor  bubble  of  "critical"  radius  and  pressure  will  be  in 
static  equilibrium  if 
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(2-1) 


2o  _  _2£ 

"  P^-P,  ■  AP 

where  o  is  the  surface  tension  of  the  liquid,  and  Pt  is  the 
pressure  of  the  saturated  liquid.  Since  bubble  detectors  are 
usually  operated  at  temperatures  sufficiently  removed  from  the 
critical  temperature,  the  pressure  difference  between  the  bubble 
(inside)  and  the  liquid  (outside)  is^ 

AP  =  P^-P,  =  (P^-P,)|l-^j  (2-2) 

Where  p,  and  p,  refer  to  the  densities  of  the  saturated  vapor  and 
saturated  liquid  respectively,  at  the  liquid's  temperature. 

Although  the  bubble-and-liquid  system  is  mechanically 
determinate  and  in  thermodynamic  equilibrium,  it  is  not 
thermodynamically  stable.  This  condition  of  instability  implies 
that  were  the  bubble  any  smaller,  it  would  collapse.  However,  if 
it  gains  one  additional  molecule,  its  radius  will  exceed  and 
grow  spontaneously  until  the  initial  liquid  droplet  has  fully 

evaporated . 

For  the  bubble  to  achieve  and  exceed  the  critical  radius,  a 
certain  critical  energy  must  be  deposited  by  the  charged 
particle  as  it  passes  through  the  superheated  liquid  droplet.  This 
minimum  amount  of  energy,  E^,  has  been  studied  at  great  length  and 
can  be  evaluated  by  examining  the  theoretical  amounts  of  energy 
demanded  by  each  of  the  mechanisms  involved  in  bubble  formation, 
critical  energy  can  be  calculated  as  the  sum  of  various 
constituents,  with  the  formula  for  E^  being 
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+  4itr^|o-r|||  +  +  P»^*  +  W’j  +  H^visc  •"  ^h 

(2-3) 

Fundamentally,  Eg  consists  of  both  reversible  and  irreversible 
components.  The  first  three  terms  represent  thermodynamically 
reversible  processes  which  account  for  the  formation  of  the  bubble 
surface,  the  evaporation  of  the  liquid  mass  inside  the  bubble  to 
the  same  mass  of  vapor,  and  the  expansion  work  done  against  the 
pressure  of  the  liquid.  The  final  four  terms  represent 
irreversible  processes  and  account  for  the  kinetic  energy  imparted 
by  the  expanding  bubble  against  the  wall  of  the  liquid,  ;  the 
acoustic  energy  lost  by  the  generation  of  sound  waves  ("popping"), 
Wj  ;  the  energy  lost  during  the  bubble's  growth  by  the  action  of 
viscous  forces,  P/viec  !  and  the  thermal  energy  lost  during  the  time 
a  bubble  is  expanding  to  its  critical  radius,  W^,.  Calculations  of 
each  term  of  Eg. (2-3)  for  various  freon  compounds  at  different 
temperatures  predict  that  the  first  three  reversible  processes 
typically  account  for  more  than  99%  of  the  critical  energy. 

All  of  the  energy  required  to  produce  the  bubble  originates 
from  the  energy  lost  by  the  charged  particle  passing  through  the 
liquid.  No  bubble  will  form  if  the  net  energy  deposited  is  less 
than  E^.  Furthermore,  E^  must  be  deposited  within  a  short  enough 
time  and  inside  a  small  enough  volume  to  be  effective.  The 
effective  length  L  over  which  a  charged  particle  transfers  its 
energy  can  be  related  to  the  dimensions  of  the  "critical  bubble." 
If  the  stopping  power  for  the  recoil  ion  is  assumed  to  be  constant, 
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'then  'the  energy  available  for  bubble  formation  is 


(2-4) 


If  the  relationship  between  L  and  were  known,  then  predicting  a 
minimum  stopping  power  necessary  to  cause  a  bubble  would  be 
possible,  assuming  Ej^  -  E^, 

Bell^*  defined  this  relationship  by  introducing  the  parameter 


a ,  such  that 


L  =  ar,  (2-5) 

Although  the  exact  value  of  L  is  undetermined,  one  might 
intuitively  think  that  the  energy  must  be  deposited  within  an 
interval  of  one  critical  bubble's  diameter,  making  a  =  2.  Table  2— 
1  shows  that  several  researchers  have  greatly  differed  in  their 
estimation  of  the  parameter  a. 


TABLE  2-1  Predicted  Values  of  Dimensionless  Parameter  a,  where 
L  =  a  r^. 


Author 

Year 

Ref. 

a 

Norman  and  Spiegler 

1963 

[19] 

2ir 

El-Nagdy  and  Harris 

1971 

[20] 

2 

Deitrich  and  Connolly 

1973 

[21] 

27r  to  12.96 

Bell 

1974 

[22] 

6.07 

Apfel 

1988 

[23] 

2 

Harper^®  however,  takes  an  unconventional  approach  and  defines  L  in 
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terms  of  a  different  parameter.  A  bubble  which  has  been  formed  in 
a  liquid  may  grow  to  visible  dimensions  if  its  radius  exceeds  a 
critical  value  r^.  Since  the  mass  of  vapor  present  in  the  critical 
bubble  has  not  changed  since  it  was  in  the  liquid  phase,  (which  is 
when  the  charged  particle  traversed  it) ,  then  it  would  make  sense 
to  redefine  the  relationship  between  effective  ion  track  length  and 
bubble  dimension  in  terms  of  this  initial  "seed”  bubble's 
dimensions.  The  smaller  seed  bubble  has  a  radius  equal  to  a 
sphere  of  liquid  having  the  same  mass  as  that  of  the  vapor  present 
in  a  critical  bubble.  Thus,  the  seed  bubble,  which  contains  this 
mass  of  vapor  before  expansion  has  begun,  has  a  radius  of 


If  the  bubble  radius  r  exceeds  when  the  ion-supplied  energy  Ej^ 
is  expended,  then  it  will  continue  to  grow  as  further  energy  is 
received  from  the  superheated  liquid.  Therefore,  Harper  argues,  it 
seems  appropriate  to  link  the  seed  bubble  radius  with  the 
effective  track  length  L  of  the  moving  ion  by  creating  the 
dimensionless  coefficient  b,  such  that 

L  =  br^  (2-7) 

where  b  has  been  empirically  shown  to  be  4.3  from  experimental 
data.^^ 

2.3  THEORETICAL  RESPONSE  MODEL 

The  response  of  a  bubble  detector  will  vary  directly  with  the 
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amount  of  superheated  liquid,  magnitude  of  neutron  flux, 
microscopic  cross  section  for  a  particular  reaction  and  the 
efficiency  factor  for  each  constituent  in  the  liquid.  The  response 
rate  R  (bubbles/sec)  can  be  calculated  in  terms  of  these  factors  by 
the  equation 


R 


where  T 

<l>iE„) 

Vt 


F/  (E„,T) 


i 

and  j 


{T,E„)  =♦{£„)  FjUt.EJ  (2-8) 

■1  j 

is  the  superheated  liquid  droplet  temperature 

is  the  incident  neutron  energy 

is  the  neutron  flux,  n/cm^/sec 

is  the  total  volume  of  superheated  liquid,  cm^ 

is  the  atomic  density  of  the  ith  species  in  the 

liquid,  atoms /cm^ 

is  the  microscopic  cross  section  for  a  j-type 
reaction  in  the  ith  species,  at  neutron  energy 
barns 

is  the  efficiency  factor  for  a  j-type  reaction  in 
the  ith  species  for  a  particular  liquid  temperature 
T  and  the  incident  neutron  energy  E^,  dimensionless 
is  one  of  the  atomic  constituents  of  the  compound 
is  the  type  of  neutron  reaction. 


The  efficiency  factor  F(E„)  refers  to  the  fraction  of  knock-on 
ions  that  are  capable  of  depositing  the  critical  amount  of  energy 
(thus  forming  a  bubble) ,  considering  their  initial  formation  energy 
spectra  and  stopping  powers  in  the  superheated  liquid  of  the 
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detector.  Because  critical  energy  depends  upon  the  superheated 
liquid's  temperature,  the  efficiency  factor  will  also  be  a  factor 
of  temperature. 

If  the  neutron  field  is  polyenergetic  —  such  as  that  found  in 
spontaneous  fission  sources  like  californium  and  nuclear  warheads  - 
-  then  the  total  response  rate  at  a  given  temperature  is  found  by 
integrating  over  all  energies 

R(T)  =  FjUT,E^)  dE  (2-9) 

E  -7 
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CHAPTER  III 


SEARCH  FOR  ALTERHATIVE  DROPLET  MATERIALS 


The  theoretical  response  model  previously  discussed  is  used  to 
identify  possible  candidate  materials  which  could  provide  a 
consistent  response  (bubbles)  with  changing  temperatures.  The 
materials  are  classified  as  exceptional  or  poor  candidates  by 
assigning  a  Figure  of  Merit  (FOM) .  For  the  three  best  candidates, 
a  response  is  calculated.  A  protocol  for  determining  constant 
responses  is  established.  The  "flat"  regions  of  response  provided 
by  the  candidate  materials  are  graphically  displayed. 


3.1  CANDIDATES  EVALUATED 

Current  superheated  liquid  neutron  detectors  are  inherently 
temperature  dependent.  Although  both  Apfel  Enterprises  and  BTI 
claim  to  have  overcome  this  problem  through  computer  algorithms  and 
pressure  compensating  devices,  the  fundamental  principle  that  the 
critical  energy  changes  with  temperature  remains  unchanged.  The 
critical  energy  is  directly  linked  to  the  themnodynamics  of  the 
bubble  formation  process  and  must  be  temperature  dependent.  The 
methods  used  by  Apfel  Enterprises  and  BTI  could  work;  however,  it 
would  be  best  to  attack  the  temperature  dependence  problem  at  its 
root  and  identify  a  possible  superheated  liquid  droplet  material 
that  would  naturally  provide  a  "flat"  detector  response  over  a 
reasonable  temperature  range.  Consequently,  a  study  was  conducted 
to  find  suitable  liquids  other  than  those  currently  employed  that 
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FIGURE  3-1 


Critical  Energy  and  bubble  radius  for  propylene  as 
a  function  of  temperature  at  1  atmosphere. 


The  evaluation  process  was  conducted  in  several  steps.  First, 
many  different  refrigerants  were  surveyed  by  studying  the 
thermodynamic  properties  and  calculating  the  critical  energy  as 
a  function  of  temperature.  Refrigerants  were  chosen  as  the 
possible  candidate  materials  because  of  their  favorable 
thermophysical  properties.  A  combination  of  low  boiling 
temperatures  at  atmospheric  pressure,  high  heats  of  vaporization, 
and  low  surface  tensions  make  refrigerants  ideal  materials  to  be 
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used  as  the  superheated  liquid  candidates  for  nucleation.  The 
critical  energy  calculations  were  performed  by  evaluating  only  the 
reversible  work  terms  of  Eg.  (2“3) .  The  modified  equation  for  Ect 

was  calculated  using  thermodynamic  properties  obtained  from  Dupont 
Laboratories,^®  Stanford,'®  and  the  National  Institute  of  Standards 
and  Technology  (NIST).'’  Thermophysical  properties  were  generally 
avalla.ble  for  temperatures  in  excess  of  310®K  but  not  for  all 


CRITICAL  BUBBLE  DATA 


iSHSBBi 


FIGURE  3-2  Critical  energy  and  bubble  radius  for  propane  as  a 

function  of  temperature  at  1  atmosphere. 
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candidates.  Theoretical  calculations,  therefore,  were  carried  out 
using  the  fullest  range  of  available  data.  Figures  3-1,  3-2,  and 
3-3  show  the  critical  energy  and  critical  radius  as  a  function 
of  temperature  for  propylene,  propane,  and  HFC-134a  respectively. 
Table  3-1  lists  the  thermophysical  properties  and  energy 
requirements  for  the  three  candidates  at  selected  temperatures.  It 
is  important  to  note  that  HFC-134a  has  a  much  higher  critical 
energy  and  boiling  temperature  than  both  propylene  and  propane. 
Data  for  some  other  compounds  studied  are  graphically  presented  in 
Appendix  A. 

The  critical  energy  needed  to  form  a  bubble  is  less  dependent 
upon  temperature  near  the  material's  critical  temperature.  Thus, 
having  a  critical  temperature  close  to  the  normal  operating 
temperatures  of  bubble  detectors  is  a  favorable  characteristic  for 
possible  superheated  liquid  droplet  candidates.  However,  by 
operating  so  close  to  the  critical  temperature,  the  issue  of 
homogeneous  nucleation,  the  spontaneous  vaporization  of  droplets 
without  neutron  interaction,  became  a  definite  consideration. 
Therefore,  the  theoretical  temperature  at  which  spontaneous 
nucleation  occurs,  the  "foam  limit,"  was  calculated  at  one 
atmosphere  pressure  using  Spiegler's  method.^®  Based  upon  the  van 
der  Waals  equation,  Spiegler's  unique  derivation  uses  only  the 
critical  constants  of  the  liquid  to  determine  the  foam  limit.  The 
guideline,^®  that  for  small  volumes  homogenous  nucleation  will  very 
rarely  (once  per  million  years)  occur  at  temperatures  more  than  3°C 
below  the  foam  limit,  was  used  when  evaluating  the  candidate 
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materials. 

After  the  critical  energies  and  foam  limit  temperatures  had 
been  calculated,  the  materials  were  classified  as  possible 
candidates  by  assigning  each  one  a  non-dimensional  Figure  of  Merit 
(FOM)  .  The  FOM  was  calculated  to  quantify  the  change  in  a 
candidate's  theoretical  for  a  given  change  in  temperature, 
within  the  normal  operating  range  of  the  detectors.  The  variation 
in  Ec  was  calculated  over  a  range  from  about  290°K  up  to  about 
310°K,  or  to  the  highest  temperature  for  which  properties  were 
available.  The  lowest  gradient  thus  calculated,  that  of  Freon— 23, 
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TABLE  3-1.  Thermophysical  Properties  and  Energy  Requirements 

for  Alternate  Bubble  Dosimeter  Compounds  at  27®C. 


Compound  Name 

Propylene 

Propane 

HFC-134a 

Chemical  Formula 

CsHe 

C3H3 

CH2FCF3 

Boiling  Point 

“C 

-47.7 

-42.07 

-26.50 

Temperature 

®c 

27 

27 

27 

Vapor  Density  (p^) 

kg/m^ 

25.5 

21.7 

32.6 

Liquid  Density  (p^) 

kg/m^ 

503.3 

489.2 

1192.9 

Vapor  Pressure 

atm 

11.9 

9.8 

6.9 

Surface  Tension  (a) 

dyne /cm 

6.4 

6.6 

7.7 

Critical  Radius 

10 -^cm 

1.22 

1.55 

2.67 

Surface  Formation 
Energy 

kev 

0.54 

0.75 

2.54 

Vaporization  Energy 

kev 

0.40 

0.70 

3.02 

Expansion  Energy 

K 

keV 

0.00 

0.01 

0.05 

Total  Critical 

Energy 

-h 

keV 

0.94 

1.46 

5.61 

was  assigned  a  FOM  of  1.00  and  all  others  were  compared  directly, 
resulting  in  some  cases  having  significantly  larger  FOM's.  It  is 
possible  to  calculate  the  theoretical  value  of  E<.  for  a  temperature 
above  the  foam  limit  and  ascertain  its  rate  of  change  with  respect 
to  temperature,  as  is  evident  for  some  cases  in  Table  3-2.  It  was 
only  after  comparing  both  the  FOM  and  foam  limit  temperatures 
however,  that  a  decision  could  be  made  concerning  which  materials 
made  the  best  candidates.  The  results  of  the  FOM  comparisons  may 
be  found  in  Table  3-2.  From  the  table  it  can  be  seen  that  although 
a  candidate  might  have  a  favorable  FOM,  the  foam  limit  temperature 
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is  often  well  below  the  temperature  region  of  interest. 
Furthermore,  the  great  disparities  in  FOM  for  the  different 
candidates  should  be  noticed.  After  evaluating  the  data,  five 
candidates  stood  out  as  potentially  offering  the  best  chance  of 
achieving  a  non-temperature  dependent  response  for  the  normal 
operating  region.  These  candidates  included:  Freon  22,  Freon  500, 
HFC-134a,  Propylene,  and  Propane. 

TABLE  3-2.  Figure  of  Merit  for  the  Candidate  Alternate 


Compounds 


Candidate 

AEe 

(keV) 

> 

O 

AEc/AT 

FOM 

Foam 

Limit 

(*K) 

Freon- 12 

19.74 

290-310 

0.99 

7534 

326 

Freon- 13 

0.001364 

290-295 

0.000273 

2.08 

256 

Freon- 2 2 

5.383 

290-310 

0.27 

2055 

312 

Freon-23 

0.000655 

290-295 

0.000131 

1.00 

253 

Freon-114 

3693 

290-310 

185 

>  10^ 

355 

Freon-115 

1.156 

290-300 

0.12 

882 

299 

Freon-C318 

150.9 

290-310 

7.55 

>  10^ 

330 

Freon-500 

7.028 

290-310 

0.35 

2682 

321 

Freon-502 

1.453 

290-310 

0.07 

>  10^ 

301 

HCFC-124 

226.5 

289-311 

10.2 

>  10^ 

335 

HFC-125 

0.740 

289-311 

0.03 

257 

287 

HFC-134a 

19.19 

289-311 

0.87 

6659 

317 

Butane 

1480 

290-310 

74 

o 

H 

A 

359 

Ethane 

0.003 

290-300 

0.0003 

2.29 

259 

Isobutane 

115.0 

290-310 

5.75 

>  10^ 

347 

Propylene 

2.286 

290-310 

0.11 

873 

309 

Propane 

3.401 

290-310 

0.17 

1298 

313 
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3.2  RESPONSE  CALCULATIONS  USING  THEORETICAL  MODEL 


Of  the  five  candidates  listed  above,  propylene  had  the 
smallest  FOM  and  was  designated  the  candidate  with  the  most 
potential  for  achieving  a  "flat"  response  over  the  normal  operating 
region.  In  order  to  establish  the  theoretical  response  for  a 
device  using  propylene  as  the  droplet  material,  the  Harper  model 
was  applied  to  a  bare  Californium  source  having  a  median  neutron 
energy  of  1.678  MeV.  The  spectrum  for  the  Californium  source  was 
broken  into  44  groups,  or  bins,  of  energy  based  on  the  standard 
groups  recommended  by  NIST.”  Individual  mean  energies  for  the 
groups  ranged  from  0.025  MeV  to  17  MeV.  The  individual  groups' 
relative  strengths  were  also  obtained  from  reference  [29] .  Figure 
3-4  illustrates  the  bare  Californium  spectrum.  The  Californium 
spectrum  was  chosen  for  the  calculations  because  it  is  a 
spontaneous  fission  neutron  source  whose  neutron  energy  spectrum  is 
well  characterized  and  useful  for  model  validation. 

The  relative  contribution  of  each  group's  flux  to  the  overall 
response,  or  group  fractional  response  (GFR) ,  was  calculated  by 
dividing  the  group's  differential  source  strength  by  the  total 
strength  given  in  reference  [29].  Dividing  the  spectra  into 
individual  bins  allows  the  conversion  of  the  energy  integral  of  Eq. 
(3-9)  into  the  summation 

/  71=44 

<!)(£■„)  dE„  = 

^  n=l 

Microscopic  elastic  scattering  cross  section  data  was  obtained 
from  Lawrence  Livermore  National  Laboratory  (LLNL)^°  and  the  BNL-325 
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BARE  Cf-252  SPECTRUM 
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FIGURE  3-4. 


Neutron  spectrum  from  bare  spontaneous  fission 
source.  The  spectrum  has  an  average  energy  of  1.678 
MeV. 


publication^^  for  each  element  of  propylene.  An  average  cross 
section  was  calculated  for  each  bin,  such  that  the  area  of  the 
rectangle  as  high  as  the  bin  and  as  wide  as  the  energy  group  was 
exactly  equal  to  the  total  area  bounded  by  the  continuous  cross 
section  versus  energy  curve,  evaluated  over  the  energy  of  that 
particular  group.  This  averaging  procedure  accounted  for  any 
possible  resonances  which  occurred  within  the  bin.  Figure  3-5 
shows  the  elastic  cross  sections  versus  neutron  energy  for  common 
refrigerant  elements  such  as  carbon,  hydrogen,  and  fluorine. 
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SCATTERING  CROSS  SECTIONS 


Neutron  Energy  (MeV) 


FIGURE  3-5.  Neutron  elastic  scattering  cross  sections  for  H,  C, 

and  F. 


Values  for  other  microscopic  cross  sections,  such  as  inelastic 
scattering  and  absorption  reactions,  were  negligible  and  were  not 
used  in  the  response  calculations. 

The  atom  density  (atoms/cm^)  used  in  the  response 

calculations  was  computed  using  the  standard  formula 

^  (3-3) 

M 

where  Nj^  is  Avogadro's  number,  p  is  the  superheated  liquid's 
density,  is  the  isotopic  abundance  of  the  particular  isotope  i, 
and  M  is  the  liquid  compound's  molecular  weight.  This  number. 
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multiplied  by  the  number  of  atoms  of  the  ith  species  per  molecule, 
gives  the  isotope's  individual  atom  density. 

The  efficiency  factor,  F(E„,T) ,  was  evaluated  for  each 
possible  recoil  ion  by  considering  both  the  recoil  ion's  initial 
energy  spectra  as  well  as  its  stopping  powers  in  the  superheated 
liquid  of  the  detector.  F(E„^T)  accounts  for  the  expected  fraction 
of  knock-on  ions  capable  of  depositing  the  critical  amount  of 
energy  within  the  detector  droplets. 

Using  Harper's  assumption  that  the  amount  of  energy  made 
available,  by  a  charged  particle  transiting  a  distance  L  =  br^ 
through  a  liquid  droplet,  is  equal  to  the  critical  energy,  E^,  an 
expression  for  minimum  stopping  power  (De/dx)  can  be  written  as 


Note  that  the  right-hand  side  of  Eq.  (3-4)  depends  only  on  the 
thermophysical  properties  of  the  superheated  liquid.  Because 
stopping  powers  can  be  calculated  as  a  function  of  an  ion's  energy, 
it  is  possible  to  determine  exactly  how  much  energy  the  individual 
ions  of  the  superheated  liquid  must  have  in  order  to  exhibit  the 
minimum  required  stopping  power,  (De/dx) Figures  3-6,  3-7,  and 
3-8  are  graphical  representations  from  Ziegler's  Transport  of  Ions 
in  Matter  (TRIM-90)  code^^  of  the  stopping  powers  for  the  possible 
recoiling  ions  produced  in  propylene,  propane,  and  HFC-134a.  For 
propylene,  3.52  MeV  per  mg/cm^  is  the  minimum  stopping  power 
required  by  Eq.  (3-4)  for  a  carbon  ion  to  form  a  bubble  at  a 
temperature  of  36.9°C.  Entering  Fig.  3-6  with  a  value  of  3.52  on 
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the  dE/dx  axis,  the  corresponding  value  for  the  reguired  carbon 
energy  is  found  to  be  approximately  0.13  MeV.  As  a  result,  only 
carbon  ions  with  more  energy  than  0.13  MeV  are  capable  of 
depositing  the  necessary  critical  energy  within  the  propylene 
liquid  droplet. 


STOPPING  POWER  IN  PROPYLENE 


FIGURE  3-6. 


Total  (compound)  stopping  power  in  propylene  for  H 
and  C  ions  as  a  function  of  ionic  energy,  generated 
from  Ziegler's  Trim-90  code. 


According  to  the  Harper  model,  all  possible  neutron 
interactions  must  be  evaluated  to  determine  whether  the  required 
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STOPPING  POWER  IN  PROPANE 


FIGURE  3-7.  Total  (compound)  stopping  power  in  propane  for  H 
and  C  ions  as  a  function  of  ionic  energy,  generated 
from  Ziegler's  Trim-90  code. 

(dE/dx)  can  be  provided  by  the  recoil  ions.  Only  those  neutron 
interactions  which  occur  within  the  droplet  are  examined  since  it 
is  extremely  unlikely  that  a  recoiling  ion  from  the  gel  matrix 
material  would  result  in  a  droplet  nucleation*.  For  the  case 
considered,  elastic  scattering  was  the  dominant  interaction  and 
resulted  in  fairly  predictable  knock-on  ion  spectra.  The  energy  of 
elastically  scattered  ions.  Eg, ,  depends  on  both  the  incident 
neutron  energy,  E^.i,  and  the  center-of-mass  scattering  angle,  6^, 
such  that  the  energy  can  be  expressed  as  a  direct  function  of  the 
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STOPPING  POWER  IN  MFC- 134a 


Ion  Energy  (MeV) 


FIGURE  3-8.  Total  (compound)  stopping  power  in  HFC-134a  for  H, 
C,  and  F  ions  as  a  function  of  ionic  energy, 
generated  from  Ziegler's  Trim-90  code. 

scattering  angle  0^  by  the  following  equation 


(^+1)2 


(1-cose^) 


(3-5) 


As  a  result,  for  a  particular  the  relative  distribution  of 
recoiling  ions  follows  a  unique  curve  related  to  the  possible 
center-of-mass  scattering  angles,  6^. 

Using  LLNL's  library  of  evaluated  nuclear  data^®,  the  expected 
elastic  scattering  distributions  were  generated  for  carbon, 
hydrogen,  and  fluorine  at  all  energies  of  interest.  Figure  3-9 
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ENERGY  DISTRIBUTION 

Elastic  Scattering,  3.50  MeV  neutrons 


Carbon  ton  Energy  (MeV) 


FIGURE  3-10.  Distribution  of  carbon  ions  after  scattering  by 

3.50  MeV  neutrons.  Shaded  area  indicates  fraction 
of  recoiling  ions  which  supply  Ec 

shows  the  relative  distribution  of  carbon  ions  scattered  by  an 
incident  3.50  MeV  neutron.  Considering  the  carbon  ions,  it  can  be 
seen  that  there  is  a  relatively  high  probability  that  6^  will  be 
close  to  either  0°  (slight  glancing  collision,  almost  a  complete 
miss)  or  180°  (head-on  collision,  complete  back-scattering) .  From 
the  ion  distribution,  it  is  possible  to  calculate  the  fractions  of 
>  ions  generated  that  have  either  energies  greater  than  a  certain 

value  or  energies  which  lie  in  a  certain  energy  band.  For 
propylene,  at  36.9°C,  it  was  earlier  stated  that  the  carbon  ions 
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FIGURE  3-9.  Relative  distribution  as  a  function  of  carbon  ion 

energy  after  elastic  scattering  3.50  MeV  neutrons, 
showing  extreme  forward  and  backward  scattering. 

needed  to  have  a  minimum  energy  of  0.13  MeV  in  order  to  form  a 
bubble.  The  fraction  of  all  carbon  ions  capable  of  creating  a 
propylene  bubble  is  found  by  comparing  the  shaded  area  of  Fig.  3-10 
to  the  total  area  under  the  relative  distribution  curve  produced  by 
elastic  scattering  with  a  3.50  MeV  neutron.  This  fraction  is  the 
efficiency  factor,  F(E„,T)  ,  and  was  calculated  to  be  0.73.  An 
efficiency  factor  of  0.73  implies  that  73%  of  the  recoil  carbon  < 

ions  created  from  elastic  scattering  with  a  3.50  MeV  neutron  have 
sufficient  energy  to  form  a  bubble.  Figure  3-11  also  demonstrates 
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the  computation  of  F(E„,T}  for  hydrogen  ions  elastically  scattered 
by  a  2.0  MeV  neutron  in  propylene  at  36.9®C.  Unlike  the  carbon 
ion,  the  hydrogen  ion  has  an  upper  bound  for  the  ion  energies  which 
can  produce  nucleation.  The  upper  energy  boundary  results  from  the 
unigue  shape  of  the  hydrogen  ion  stopping  power  curve.  From  Fig. 
3-6,  it  should  be  observed  that  hydrogen  ion  energies  greater  than 
0.1  MeV  begin  to  provide  a  decreasing  (dE/dx)  with  increasing  ion 
energy,  after  the  Bragg  peak  energy  of  about  0.1  MeV  (100  keV) .  As 
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Elastic  Scattering,  2.0  MeV  Neutrons 


Hydrogen  Ion  Energy  (MeV) 


FIGURE  3-11.  Distribution  of  hydrogen  ions  after  scattering  by 

2.0  MeV  neutrons.  Shaded  area  indicates  fraction 
of  recoiling  protons  which  supply  E^. 

demonstrated  in  Fig.  3-6,  some  higher  energy  hydrogen  ions  are  not 
capable  of  supplying  a  (dE/dx)  which  lower  energy  ions  were  able  to 


do. 

In  determining  the  propylene  response  for  a  bare  Californium 
spectrum,  efficiency  factors  were  calculated  for  elastic  scattering 
with  neutrons  from  each  of  the  44  different  energy  bins,  having 
energies  ranging  from  0.0253  eV  to  14  MeV.  Furthermore,  F (E„,T)va.s 

1 

evaluated  at  temperatures  ranging  from  -13.2®C  to  46.9®C.  Figure 
3-12  shows  how  the  F(E„,T)  for  carbon  ions  varies  for  a  range  of 
temperatures  and  neutron  energies  within  propylene.  As  seen  in  the 
plot,  the  higher  the  incident  neutron's  energy  or  the  higher  the 
liquid's  temperature,  the  closer  the  efficiency  factor  will  be  to 
unity . 

This  is  the  normal  trend  for  most  ions.  However,  as  shown  in 
Fig.  3-13,  the  F(E„,T)  for  hydrogen  ions  behaves  quite  differently. 
Because  hydrogen  ions  are  not  capable  of  providing  a  large  (dE/dx) , 
the  stopping  power  curve  is  very  unique  when  compared  with  other 
ions.  As  the  ion  energy  increases  above  1  MeV,  the  (dE/dx) 
provided  by  the  ion  becomes  too  small  to  form  a  bubble.  As  a 
result,  the  efficiency  factor  decreases  with  increasing  neutron 
energy.  However,  because  the  required  (dE/dx)  decreases  with 
increasing  temperature,  the  efficiency  factors  always  increase  with 
increasing  temperature.  When  F (E, „^T)=1 ,  then  every  recoil  ion  is 
capable  of  producing  a  bubble.  Consequently,  the  detector  response 
will  be  "flat"  when  the  total  efficiency  factor  for  all 
recoil  ions  achieves  unity,  < 

A  close  evaluation  of  the  stopping  power  plots  for  carbon  and 
hydrogen  ions  in  propylene  reveals  that  there  are  threshold 
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EFFICIENCY  FACTOR 

Carbon  Ions  in  Propylene 
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FIGURE  3-12.  Efficiency  factor  F(E„,T)  for  carbon  ions  in 
propylene  as  a  function  of  temperature. 

stopping  powers  below  which  all  recoil  ions  of  a  particular  type 
can  satisfy  the  critical  energy  requirements.  As  a  result,  if  the 
required  (dE/dx)  is  less  than  the  threshold  for  a  particular  ion, 
recoil  ions  of  all  energies  will  have  more  stopping  power  than  the 
minimum  required,  and  thus  be  capable  of  supplying  the  necessary 
critical  energy,  E^.  For  the  case  of  propylene,  as  seen  in  Fig.  3- 
6,  if  (dE/dx)  is  less  than  2  MeV  per  mg/cm%  then  recoiling  carbon 
ions  of  all  energies  are  capable  of  providing  the  critical  energy. 
As  a  result,  F(E„,T)  will  equal  one  for  carbon  ions  of  all 
energies.  The  propylene  response  will  continue  to  be  constant 
until  (dE/dx) decreases  below  the  upper  threshold  (largest  dE/dx 
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EFFICIENCY  FACTOR 

Hydrogen  Ions  in  Propylene 
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FIGURE  3-13.  Efficiency  factor  F(E„^T)  for  hydrogen  ions  in 
propylene  as  a  function  of  temperature. 


provided  by  ion)  for  hydrogen  ions.  Previously,  hydrogen  ions  of 
all  energies  had  been  incapable  of  supplying  the  required  (dE/dx) . 
Once  the  upper  threshold  is  passed,  hydrogen  ions  will  begin  to 
contribute  to  the  response,  resulting  in  a  total  response  that  is 
no  longer  flat.  The  propylene  response  cannot  become  constant 
again  until  the  required  (dE/dx)  becomes  less  than  the  lower 
hydrogen  threshold  (dE/dx  below  which  all  resulting  ions  are 
capable  of  depositing  E^)  .  However,  when  this  occurs,  the 
temperature  would  have  increased  so  much  that  the  material  would 
have  already  reached  its  foam  limit  of  35.7°C.  Figure  3-14  shows 
the  propylene  response  for  the  polyenergetic  neutrons  of  bare 
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PROPYLENE  RESPONSE 

Bare  Californium 


FIGURE  3-14.  Theoretical  model's  predicted  response  for 
propylene  (CjHg)  to  bare  Californium  neutron 
spectrum  at  1  atmosphere. 

Californium.  The  response  has  several  important  characteristics 
that  should  be  noted.  First,  there  are  no  hydrogen  effects  until 
24 ®C  because  at  lower  temperatures,  (dE/dx)^^^  is  too  high  and 
cannot  be  supplied  by  any  energy  hydrogen  ion.  It  is  not  until 
temperatures  greater  than  24 ®C  are  reached  that  is  low  enough 
for  (dE/dx)„i„  to  be  satisfied  by  hydrogen  ions.  Second,  one  can 

■m 

see  the  F(E„,T)  curves  echoed  in  the  individual  ion's  responses  as 
^  expected  by  Eg.  (2-9) .  For  carbon,  all  possible  ions  have  a 

F(E„fT)=l  at  a  temperature  of  16.8°C. 

The  theoretical  responses  for  2.06  and  2.1  MeV  neutrons. 
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respectively,  were  calculated  for  propane  and  HFC-134a.  As  shown 
in  Figs.  3-15  and  3-16  respectively,  plateau  regions  of  relatively 
constant  response  also  exist  in  unique  temperature  bands  for  these 
other  alternate  materials.  Although  the  response  was  only 
calculated  for  a  2.1  MeV  neutron  source  and  not  for  every  energy 

1 

bin,  it  still  gives  a  good  indication  of  the  trends  of  the  complete 
response.  The  constant  temperature  plateau  should  span  the  same 
temperature  band,  only  the  response  before  and  after  the  plateau 
should  change.  The  boundaries  of  the  response  plateau  are 
determined  more  by  (dE/dx)  than  the  incident  neutron  energy.  A 
2.1  MeV  neutron  has  sufficient  energy  so  that  at  least  a  portion  of 
the  knock-on  hydrogen  ions  created  will  have  enough  energy  to  form 
a  bubble  if  (dE/dx)  is  less  than  the  upper  threshold  for 
hydrogen.  Because  the  response  only  considers  2.1  MeV  neutrons, 
the  response  before  and  after  the  plateau  does  not  account  for  the 
combined  responses  of  the  other  43  energy  bins  as  dictated  by  Eq. 

(2-9) . 

3 . 3  CONCLUSIONS 

The  extent  of  the  regions  of  constant  response  is  directly 
related  to  both  the  required  stopping  power,  (dE/dx) of  each 
material  and  the  stopping  power  thresholds  of  the  individual 
components  of  the  materials.  This  is  very  clear  when  considering 
propane  and  propylene.  For  these  simple  hydrocarbons,  the  plateau  ^ 

begins  once  the  required  (dE/dx)  falls  below  the  minimum  lower 
threshold  for  the  recoil  carbon  ions.  At  this  point,  the  upper 
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PROPANE  RESPONSE 

2.06  MeV  Neutron 


FIGURE  3-15.  Theoretical  model's  predicted  response  for  propane 

(CjHg)  to  bare  Californium  neutron  spectrum  at  1 
atmosphere . 

threshold  for  a  hydrogen  response  is  significant  when  compared  with 
the  required  stopping  power.  However,  as  the  temperature 
increases,  the  required  (dE/dx)  gradually  decreases  in  value  and 
eventually  becomes  less  than  the  upper  threshold  of  the  hydrogen 
recoil  ions.  When  this  occurs,  the  response  is  no  longer  flat 

* 

because  the  combined  response  is  a  summation  of  both  the  flat 
>  carbon  response  and  the  increasing  hydrogen  response.  HFC-134a  has 

a  slightly  more  complicated  response.  In  order  for  the  response  to 
be  flat,  the  required  (dE/dx)  must  be  less  than  the  lower  threshold 
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FIGURE  3-16, 


Theoretical  model's  predicted  response  for  HFC-134a 
(CH2FCF3)  to  bare  Californium  neutron  spectrum  at  1 
atmosphere. 


for  both  carbon  and  fluorine.  Once  again,  the  overall  response 
will  no  longer  be  constant  once  the  required  (dE/dx)  falls  below 
the  upper  threshold  for  hydrogen  ions. 

As  a  result  of  this  alternatives  study,  a  specific  protocol 
has  been  established  in  order  to  determine  if  a  candidate  material 
will  have  a  constant  response  over  a  specific  temperature  band. 
First,  a  candidate  material's  (dE/dx) must  be  evaluated  over  the 
operating  region  of  interest.  If  (dE/dx) is  of  small  value  and 
changes  little  over  the  specific  temperature  band,  then  the 
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material  is  classified  as  a  possible  candidate.  Second,  the  foam 
limit  temperature  must  be  checked  to  see  if  it  lies  within  the 
temperature  range  of  interest.  If  (dE/dx)„i„  is  relatively  small 
and  constant,  then  the  material  is  probably  close  to  its  critical 
temperature,  and  the  foam  limit  temperature  is  of  concern.  Third, 
the  stopping  power  curves  for  the  component  ions  within  the 
candidate  material  must  be  evaluated  for  any  lower  thresholds  below 
which  all  ions  of  a  certain  kind  can  deposit  The  response  for 
an  individual  constituent  will  be  constant  if  (dE/dx)„i„  is  less 
than  the  lower  threshold  for  that  constituent.  If  the  responses 
for  each  constituent  are  either  constant  or  zero,  then  the  combined 
response  will  also  be  constant. 

Although  the  temperature  band  for  a  constant  response  is 
usually  limited  to  about  10®C,  devices  theoretically  can  be 
manufactured  to  meet  selected  temperature  ranges  by  using  different 
alternate  materials.  If  one  knows  approximately  what  the 
anticipated  operating  temperature  band  will  be,  one  could  choose 
the  appropriate  device  so  that  a  detector  with  a  flat  response  is 
used.  Figure  3-17  demonstrates  this  concept.  In  this  plot,  the 
relative  responses  of  four  different  devices  are  superimposed  upon 
each  other.  BTI's  commercially  available  uncompensated  BD-IOOR 
clearly  demonstrates  a  non-linear  response  versus  temperature. 
However,  the  three  alternative  material  devices  considered  in  this 
report  each  have  a  unique  region  of  constant  response.  Clearly, 
there  is  a  particular  material  for  a  specific  temperature  band. 
The  flat  region  which  has  a  change  in  response  of  plus  or  minus 
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DETECTOR  RELATIVE  RESPONSES 


FIGURE  3-17.  The  predicted  relative  response  as  a  function  of 
temperature  for  a  propylene,  propane,  and  HFC-134a 
bubble  detector,  compared  to  a  BD-IOOR  detector. 


five  percent  for  propane,  propylene,  and  HFC-134a  is  approximately 
7.5°C  to  IS.l^C,  15.9°C  to  24.0°C,  and  27.0°C  to  44.2°C,  as 
illustrated  in  Fig.  3-18.  When  evaluating  Fig.  3-18,  it  should  be 
noted  that  the  bars  for  propane  and  HFC-134a  only  consider  a  2.1 
MeV  neutron  response.  Because  the  regions  were  determined  by  a 
change  in  response  of  plus  or  minus  five  percent,  they  are  defined 
beyond  the  response  plateau.  Consequently,  the  actual  "flat” 
regions  for  propane  and  HFC-134a  will  be  smaller  when  considering 
the  response  for  a  poly-energetic  source. 
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HFC-I34a 
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FIGURE  3-18. 


CHAPTER  IV 


VERIFICATION  OF  THE  THEORETICAL  MODEL 


The  theoretical  model  is  verified  through  experimentation 
conducted  at  the  National  Institute  of  Standards  and  Technology 
(NIST) .  The  experimental  set-up  is  presented.  Test  results  are 
presented  and  discussed. 


4 . 1  PROPOSED  EXPERIMENT 

The  theoretical  response  predications  presented  in  Chapter  III 
depend  heavily  upon  the  validity  of  the  Harper  model.  Harper  has 
provided  evidence  as  to  the  validity  of  his  theory  through  prior 
experimentation®.  During  this  research  effort,  the  opportunity 
arose  to  gather  additional  verification,  while  testing  the 
temperature  response  characteristics  of  the  Apfel  bubble 
dosimeters.  In  the  proposed  experiment,  Freon-12  detectors 
produced  by  Apfel  Enterprises  are  irradiated  at  various 
temperatures  using  a  pure  thermal  neutron  beam  (0.0253  Ev)  from  the 
Research  Reactor  at  NIST.  The  detector  responses  at  the  different 
temperatures  are  then  compared  to  the  predictions  of  the 
theoretical  model  in  order  to  test  the  model's  validity. 

4 . 2  THEORY 

Because  the  efficiency  factors  for  elastic  scattering 
reactions  by  thermal  neutrons  in  Freon-12  are  zero,  there  is  no 
expected  response  for  any  temperature  from  this  type  of  reaction. 
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FIGURE  4-1.  Estimated  thermal  shutoff  of  Freon-12  based  neutron 

bubble  detector  based  upon  theoretically  predicted 
critical  energy, 

However,  there  is  a  non-negligible  (0.4  barn)  thermal  neutron  cross 
section  for  the  Cl^®(n,p)S^®  reaction.  This  particular  reaction 
liberates  615  keV  of  energy,  divided  between  the  heavier  sulfur  ion 
(17  keV)  and  the  lighter  proton  (598  keV) .  Although  it  contains 
more  energy,  the  proton's  stopping  power  in  Freon-12  is  only  9%  of 
sulfur's  (0.33  keV  per  10  ®  cm  as  compared  to  3.68  keV  per  10'®  cm) 
and  thus,  is  unable  to  deposit  more  than  about  1  keV  within  the 
effective  ionic  energy  transfer  length  L.  Because  the  thermal 
neutron  absorption  reaction  with  chlorine— 35  is  the  only  possible 
source  of  energy,  the  total  amount  of  energy  made  available  to  the 
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superheated  liquid  can  be  calculated  as  the  sum  of  the  energies 
deposited  by  both  ions,  approximately  18  keV.  Using  Eg.  (3-1) , 

*  47tr|  (3-1) 

the  critical  energy  required  to  form  a  bubble  at  different 
temperatures  can  be  calculated.  Table  4-1  provides  the  critical 
bubble  data  for  Freon-12  at  some  temperatures  of  interest. 


TABLE  4-1.  Critical  Bubble 

Temperatures . 

Data  for 

Freon- 

•12  at 

Selected 

Freon-12 

(CCl^Fj) 

Temperature  “C 

18.9 

19.1 

20.0 

Vapor  Density 

kg/m^ 

31.5 

31.6 

32.4 

Liquid  Density  (p,) 

kg/m^ 

1334 

1333 

1330 

Surface  Tension  (a) 

dyne /cm 

9.9 

9.8 

9.7 

Critical  Radius 

lO'®  cm 

4.51 

4.47 

4.28 

Total  Critical  Energy 

keV 

18.2 

17.8 

16.1 

Figure  4-1  illustrates  the  critical  bubble  data  found  in  Table 
4-1  and  indicates  that  if  18  keV  of  energy  were  deposited,  it  would 
be  sufficient  to  cause  vaporization  of  Freon-12  droplets  at  a 
temperature  greater  than  or  equal  to  19 °C.  Thus,  by  irradiating 
Freon-12  detectors  at  temperatures  below  and  above  19°C,  it  is 
possible  to  validate  the  theoretical  19  °C  shutoff  /  turn-on 
temperature . 
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4 . 3  EXPERIMENTAL  SET-UP 

The  detectors  used  in  the  NIST  experiments  utilize  superheated 
liquid  Freon-12  droplets  uniformly  distributed  throughout  a  semi- 
viscous  gel  matrix.  The  detector  system,  manufactured  by  Apfel 
Enterprises,  uses  a  personal  computer  driven  device  that  holds  a 
glass  vial  containing  approximately  30,000  of  the  freon  liquid 
droplets.  Each  droplet  measures  about  lOO/xm  in  diameter  and  forms 
a  bubble  of  about  650/im  diameter  when  activated  by  sufficient 
energy.  When  the  individual  droplets  vaporize,  they  produce  a 
pressure  pulse  lasting  15—30  msec  with  a  frequency  spectrum 
primarily  in  the  IkHz-lOkHz  range.  A  piezoelectric  acoustic 
transducer  is  coupled  to  the  glass  vial  in  order  to  detect  each 
pulse.  The  output  is  electronically  filtered  to  discriminate 
authentic  neutron  events  from  any  background  noise  and  then  is 
downloaded  on  to  a  personal  computer  via  the  Apfel  detector. 

Detector  responses  were  measured  at  the  research  reactor 
located  in  NIST.  The  thermal  neutron  irradiations  were  performed 
using  the  lO-inch  filtered  pure  thermal  beam  of  the  research 
reactor.  The  filtered  thermal  beam  had  a  cadmium  ratio  of  10“  and 
provided  1.2  x  10“  neutrons  per  cm^-sec  (±3%)  for  all  irradiations. 
The  irradiations  were  performed  with  the  devices  held  at  various 
temperatures  maintained  by  a  Fisher  Scientific  low  temperature 
incubator.  Figure  4-2  depicts  the  experimental  setup. 

4.4  EXPERIMENTAL  RESULTS 

Six  of  Apfel 's  Freon-12  based  SD-100  devices  were  irradiated 


61 


FIGURE  4-2.  Schematic  of  the  experimental  setup  at  NIST.  The 

incubator  was  located  approximately  five  feet  from 
the  entrance  of  the  thermal  neutron  beam. 

using  the  NIST  thermal  beam.  The  irradiations  lasted  from  ten  to 
twenty  minutes,  and  were  conducted  at  five  different  temperatures 
ranging  from  14 ®C  to  26 ®C.  The  measured  response  rates 
(Bubbles/sec)  are  shown  in  Fig.  4-3.  It  is  important  to  note  the 
distinct  lack  of  response  below  19 °C  and  the  noticeable  increase  in 
detector  sensitivity  above  this  shutoff  point.  Although  the 
detectors'  responses  were  significantly  different  from  one  another 
due  to  differing  sensitivities,  each  detector  demonstrated  a 
definite  shutoff  at  approximately  19®C. 
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SDD-100  (FREON-12)  RESPONSE 

NIST  12/92  -  MJ.  Harper  and  J.C.  Rich 


FIGURE  4-3.  Experimental  results  for  Freon-12  (SDD-100) 
detectors  exposed  to  thermal  neutron  beam. 
Detectors  effectively  shutoff  at  19®C. 

4.5  EXPERIMENTAL  CONCLUSIONS 

The  thermal  irradiation  experiments  conducted  with  the  Freon- 
12  based  detectors  proved  that  the  model's  prediction  of  the 
minimum  critical  energy  required  to  nucleate  a  bubble  is  valid. 
From  Eq.  (3-1)  the  predicted  shutoff  temperature  for  thermal 
response  by  Freon-12  superheated  liquid  detectors  was  19 °C,  exactly 
the  same  as  that  observed  during  the  tests.  Although  Harper  has 
done  so  in  previous  experimentation,  this  set  of  tests  could  do 
nothing  to  validate  the  second  criterion  for  nucleation,  that  dE/dx 
be  such  that  (dE/dx)  x  L  ^  E^.  It  does  not  matter  whether  one 
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assumes  L  =  2  x  or  L  -  4.3  x  r^,  since  the  shutoff  temperature 
for  this  criterion  to  be  satisfied  is  lower  than  19  •’C.  At  this 
point  the  detectors  have  ceased  to  function  because  the  energy 
supplied  by  the  charged  particles,  is  less  than  the  minimum 
critical  energy  required, 
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CHAPTER  V 


EXPERIMENTAL  EVALUATIONS  OF  THE  BUBBLE  DETECTOR 

Several  series  of  experiments  were  conducted  at  the  U.S.  Naval 
Academy's  neutron  generator  facility  to  evaluate  the  response 
linearity,  repeatability,  bubble  growth  rate,  and  effects  of  aging 
on  the  neutron  bubble  detector.  The  experimental  set-ups  are 
provided  and  the  results  are  discussed.  Conclusions  are  drawn  as 
to  the  suitability  of  the  bubble  detectors  as  on-site  treaty 
verification  tools. 

5 . 1  REPEATABILITY  STUDY 

There  are  currently  two  different  manufacturers  of  bubble 
detector  devices,  Apfel  Enterprises  and  BTI  Industries  (BTI) . 
Although  both  types  of  devices  perform  on  the  same  theoretical 
principles,  their  characteristics  and  modes  of  operation  are  quite 
different.  One  of  the  greatest  differences  between  the  devices  is 
the  inability  of  the  Apfel  device  to  be  reused.  The  viscous  gel 
matrix  which  holds  the  initial  liquid  droplets  is  not  rigid  and 
consequently,  is  unable  to  keep  the  vapor  bubbles  from  joining 
together  once  they  have  formed.  As  a  result,  it  is  impossible  to 
return  the  vapor  back  into  the  original  liquid  droplets.  The  BTI 
device  however,  utilizes  a  rigid  matrix  which  keeps  the  newly 
formed  vapor  bubbles  in  place.  After  the  bubbles  have  formed  and 
the  detector  has  been  read,  it  is  relatively  easy  to  re-condense 
the  bubbles  back  into  their  original  liquid  form. 
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since  the  BTI  detector  is  reusable,  this  would  result  in  lower 
life  cycle  cost.  However,  for  the  reusable  detector  to  be  of  any 
value,  the  device's  sensitivity  must  remain  constant  after  repeated 
use.  In  the  following  experiment,  several  of  BTI's  temperature 
compensated  BD-IOOR  devices  were  repeatedly  irradiated  in  order  to 
determine  the  effect  of  reuse  on  device  sensitivity. 


EXPERIMENTAL  SET-UP 


Note:  Some  tests  used  fewer  or  more  than  the  ten  detector  positions  shown  here. 


FIGURE  5-1.  Representative  schematic  of  the  experimental  setup 

for  tests  conducted  at  USNA  neutron  generator 
facility. 
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The  BTI  BD-IOOR  bubble  detector  responses  were  measured  in  the 
Nucleonics  Lab  located  in  Rickover  Hall.  Neutron  irradiations  were 
performed  using  a  bare  Californium  source  (Cf*”)  which  provided  an 
emission  rate  of  6.136x10*  neutrons  per  second.  Fourteen  BD-IOOR 
detectors  were  attached  and  evenly  spaced  on  a  circular  aluminum 
ring  having  a  diameter  of  41.3  cm.  As  seen  in  Fig.  5-1,  the 
Californium  was  placed  in  the  exact  center  of  the  ring  so  that  each 
device  would  receive  the  same  neutron  exposure  from  the  isotropic 
source  during  each  successive  test.  The  neutron  irradiations 
lasted  exactly  five  minutes  for  each  run.  An  irradiation  time  of 
five  minutes  was  chosen  because  calculations  indicated  that  this 
exposure  time  would  produce  s  statistically  significant  number  of 
bubbles,  and  average  out  random  fluctuations.  All  irradiations 
were  conducted  at  standard  atmospheric  pressure  and  constant  room 
temperature  (22 "C) .  The  bubbles  were  counted  using  a  BTI  BDR- 
Series  II  Bubble  Reader,  a  diagram  of  which  is  shown  in  Fig.  5-2. 
The  reader  was  set  up,  configured,  and  operated  in  accordance  with 
the  procedures  delineated  in  the  BTI  Bubble  Reader  BDR-Series  II 
User's  Manual^*,  of  June  11,  1990,  as  revised  by  Rev.  3  of  July  3, 
1991. A  full  description  of  the  optical  reader  is  provided  in 
Section  5.5. 

DISCUSSION  OF  RESULTS 

In  this  study,  each  device  was  irradiated  ten  separate  times. 
Figures  5-3,  5-4,  and  5-5  provide  a  graphical  display  of  the  data 
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FIGURE  5-2.  Block  diagram  of  the  BTI  BDR-Series  II  Optical 


Reader.  The  reader  uses  two  cameras  and  a  computer 
to  count  the  bubbles. 


acquired  for  each  of  the  devices,  which  are  labeled  A  through  N. 
In  these  plots,  the  normalized  relative  response  of  each  device  is 
given  for  individual  runs.  The  normalized  response  was  obtained  by 
dividing  the  observed  number  of  bubbles  for  a  given  run  by  the 
average  number  of  bubbles  observed  over  all  exposures.  The  shape 
of  the  individual  three-dimensional  relative  response  plots  provide 
an  indication  of  how  the  bubble  counts  for  a  particular  run 
compared  to  the  average.  The  flatness  of  the  ribbons  is  indicative 
of  the  response  consistency  with  each  consecutive  run. 

Since  radiation  interactions  are  random  in  nature,  any 
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FIGURE  5-3.  Three-dimensional  display  of  the  data  obtained  in 

the  repeatability  test  for  devices  A  through  E. 
Flatness  indicates  consistency  of  response. 

observation  is  subject  to  some  degree  of  statistical  fluctuation. 
These  inherent  fluctuations  provide  a  means  to  check  the  normal 
functioning  of  any  piece  of  nuclear  counting  equipment.  If  the 
amount  of  fluctuation  associated  with  the  collected  data  is  not 
consistent  with  the  predictions  of  statistical  tests,  it  may  be 
surmised  that  an  anomaly  exists  in  the  detection  system.  Since 
nuclear  events  are  typically  described  by  a  Poisson  distribution, 
the  "Chi-squared  Test"  (x^)  can  be  used  to  determine  if  the 
experimental  data  follows  this  distribution. 

For  this  experiment,  the  two  percent  x"  test  was  applied  to 
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FIGURE  5-4.  Three-dimensional  display  of  the  data  obtained  in 

the  repeatability  test  for  devices  F  through  J. 
Flatness  indicates  consistency  of  response. 

the  data  collected.  The  test  indicated  that  the  fluctuations  in 
the  data  for  all  devices,  except  one,  fell  within  the  normal  bounds 
expected.  Device  A  failed  the  test  because  it  exhibited  abnormally 
small  fluctuations  (i.e.  its  response  did  not  exhibit  sufficient 
statistical  variation) .  However,  this  result  was  not  unexpected, 
since  statistically  one  device  in  25  is  expected  to  fail  the  two 
percent  test.” 

Counting  statistics  were  further  used  to  give  an  indication  of 
the  consistency  and  uncertainty  associated  with  the  data  obtained 
from  the  fourteen  devices.  Figure  5-6  provides  a  measurement  of 


FIGURE  5-5.  Three  dimensional  display  of  the  data  obtained  in 

the  repeatability  test  for  devices  K  through  N. 
Flatness  indicates  consistency  of  response. 

the  percent  standard  deviation  in  the  data  for  each  device.  In 
this  figure,  it  can  be  seen  that  the  percent  standard  deviation 
ranges  from  approximately  two  percent  to  7.5  percent.  Since  the 
average  number  of  bubbles  formed  for  each  run  was  365,  it  was 
expected  based  on  Poisson  statistics  that  approximately  one  third 
of  the  devices  would  have  a  percent  standard  deviation  greater  than 
5.2  percent.  As  can  be  seen  in  Fig.  5—6,  five  of  the  fourteen 
devices,  approximately  one— third,  had  percent  standard  deviations 
greater  than  5.2  percent.  It  appears  from  this  data,  therefore, 
that  the  sensitivity  of  the  devices  remains  essentially  the  same 
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for  up  to  ten  exposures,  suggesting  that  they  do,  indeed,  produce 
repeatable  results  at  these  exposure  levels. 


5.2  BUBBLE  GROWTH  STUDY 

After  a  superheated  liquid  droplet  vaporizes,  the  resulting 
bubble  formed  in  the  BD-IOOR  is  not  in  dynamic  equilibrium  with  the 
rigid  gel  matrix.  The  bubbles  continue  to  grow  after  their  initial 
formation.  Studies  conducted  with  the  early  models  of  the  BD-IOOR 
indicated  that  the  measured  number  of  bubbles  changed  depending  on 
how  long  after  irradiation  they  were  read.^'*  It  was  not  until  more 
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than  ten  hours  after  irradiation  that  the  devices  could  be  read 
with  confidence  that  the  exact  number  of  bubbles  had  stabilized  and 
was  accurate.  It  would  be  best  if  the  devices  could  be  reliably 
read  as  soon  as  possible  after  irradiation  so  that  on-site  feedback 
could  be  obtained.  Furthermore,  the  question  as  to  how  much  time 
must  elapse  before  the  devices  are  read  with  confidence  should  be 
answered . 

The  newer  temperature-compensated  BD-IOOR  devices  are  quite 
different  from  their  predecessors.  The  processing  of  the  liquid 
droplets  and  the  gel  has  been  changed  in  the  newer  devices. 
Furthermore,  the  manufacturer,  claims  that  the  bubbles  are  visible 
immediately  upon  nucleation.  Consequently,  it  was  determined  that 
the  phenomenon  of  bubble  growth  should  be  studied  in  greater  depth 
than  heretofore.  This  study  examined  bubble  growth  in  two  steps. 
The  first  step  involved  quantifying  the  actual  change  in  size  of  a 
bubble  over  time  by  measuring  the  changing  bubble  diameter.  In  the 
second  step,  the  devices  were  continuously  read  over  time  and  the 
change  in  response  (actual  number  of  bubbles  registered)  was 
evaluated. 

EXPERIMENTAL  SET-UP 

For  step  one,  the  bubble  physical  dimensional  study,  two 
different  BD-IOOR  devices  were  each  irradiated  for  exactly  ten 
seconds  using  the  Cf^®^  source  mentioned  earlier.  Each  detector  was 
exposed  to  the  same  radiation  field  and  each  radiation  was 
conducted  at  standard  atmospheric  pressure  and  constant  room 
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temperature  (22®C)  .  Any  change  in  the  size  of  the  bubble  was 
observed  by  measuring  the  actual  diameter  of  a  single  bubble  over 
time  using  the  Buehler  Omnimet  Image  Analysis  System,  at  the  Naval 
Academy's  materials  testing  laboratory. 

In  the  second  step  of  the  experiment,  three  separate  groups  of 
BD-IOOR  devices  were  irradiated  using  different  neutron  sources. 
The  first  group  of  ten  BD-IOOR  devices  was  independently  irradiated 
twice  for  thirty  seconds  using  the  Cf^®^  source.  Thirty  seconds  was 
chosen  because  approximately  sixty  bubbles  would  be  produced  in 
each  device.  This  number  of  bubbles  can  be  counted  very  accurately 
with  the  BTI  BDR-Series  II  Reader  Evaluation  System.  To  evaluate 
the  effects  of  changes  in  bubble  diameter  upon  number  of  bubbles 
counted,  the  devices  were  counted  every  30  seconds  between  one  and 
six  minutes  after  irradiation  and  every  minute  between  six  and 
twelve  minutes  after  irradiation.  The  second  group  of  ten  BD-IOOR 
devices  was  independently  irradiated  once  for  thirty  seconds  using 
a  Plutonium-Beryllium  (Pu-Be)  source.  For  this  source,  the  Pu”*  in 
the  Pu-Be  source  emits  alphas  which  react  with  the  beryllium  to 
produce  neutrons  with  an  approximate  average  energy  of  4.0  MeV. 
The  source  provided  an  emission  rate  of  2  x  10®  neutrons  per 
second.  For  this  test,  the  devices  were  read  every  30  seconds 
between  one  and  six  minutes  after  irradiation  and  every  minute 
between  six  and  twelve  minutes  after  irradiation.  Finally,  a  group 
of  four  BD-IOOR  devices  was  irradiated  for  30  seconds  using  the 
Cf^®^  source.  These  devices  were  read  periodically  over  an  1880 
minute  period  so  as  to  determine  the  effects  upon  long  term 
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coiinting.  Each  set  of  irradiations  for  the  bubble  growth  analysis 
was  conducted  under  the  same  conditions.  Furthermore,  each  group 
of  irradiations,  except  for  minor  variations  (different  source, 
different  number  of  detectors) ,  used  the  same  experimental  set-up 
as  that  outlined  in  Fig.  5-1.  The  goal  was  to  discover  whether 
increasing  bubble  sizes  resulted  in  increased  numbers  of  bubbles 
counted . 

DISCUSSION  OF  RESULTS 

The  results  of  the  first  step  are  presented  in  Fig.  5-7,  where 
the  changing  bubble  diameters  are  plotted  with  respect  to  time. 
From  Fig.  5-7,  it  appears  that  the  initial  rate  of  growth  is  quite 
large,  and  with  time,  the  rate  of  growth  begins  to  decrease. 
Although  the  curves  for  the  different  bubbles  are  not  exactly  the 
same,  it  appears  that  the  growth  of  the  two  bubbles  follows  the 
same  trend.  Differing  growth  rates  can  be  rationalized  by  various 
factors  such  as  slight  deviations  in  temperature  and  dissimilar 
liquid  droplet  densities  within  the  gel  matrix.  The  growth  is  very 
large  initially  and  then  tapers  off  with  time. 

The  initial  growth  study  indicates  that  it  takes  a 
considerable  amount  of  time,  more  than  10  days,  for  a  bubble  to 
begin  to  reach  a  state  of  equilibrium.  It  is  clear  that  the 
initial  rate  of  growth  is  considerable,  but  the  question  of  how 
much  of  an  effect  the  growth  rate  has  upon  the  reader's  ability  to 
detect  and  count  accurately  a  still-growing  bubble  must  be 
determined . 
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bubble  growth  is  rapid. 

The  second  step  of  the  bubble  growth  analysis  involved 
irradiating  different  groups  of  bubble  detectors  with  different 
sources  and  reading  them  at  various  times  after  irradiation.  The 
goal  of  the  study  was  to  determine  if  the  number  of  bubbles  counted 
by  the  optical  reader  varied  significantly  with  time  after 
irradiation,  and  also  whether  a  different  energy  neutron  source 
yielded  an  effect  upon  those  readings. 

The  first  group  of  ten  devices  was  irradiated  twice  using  the 
Cf^“  source.  The  experimental  results  of  the  first  two  runs  were 
analyzed  using  the  concept  of  the  "percent  change  from  horizontal." 
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Figure  5-8  provides  a  graphical  representation  of  this  concept, 
using  the  data  from  one  of  the  detectors.  In  this  figure,  the 
normalized  relative  response  is  plotted  with  respect  to  time  after 
irradiation.  The  response  was  normalized  by  dividing  the  number  of 
bubbles  counted  at  a  specific  time  by  the  average  number  of  bubbles 
for  all  counts.  A  linear  regression  was  performed  on  the  data  so 
that  the  best  fit  line  which  approximated  the  data  could  be 
determined.  The  percent  change  from  the  horizontal  was  next 
defined  by  multiplying  the  slope  of  the  best  fit  line  by  100 
percent.  This  calculation  provides  a  number  which  represents  the 
change  in  the  bubble  counts  with  time  in  terms  of  a  percent 
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relative  increase  or  decrease  from  the  horizontal.  Thus,  a  percent 
change  from  the  horizontal  of  zero  would  indicate  a  response  which 
does  not  change  with  the  amount  of  time  since  irradiation,  while  a 
percent  change  above  the  horizontal  of  one  would  indicate  a 
response  which  increased  by  one  percent  for  every  minute  since 
irradiation. 

Figures  5-9  and  5-10  provide  the  percent  changes  from  the 
horizontal  for  the  two  runs  of  the  devices  irradiated  by  the  Cf^®^ 
source.  The  changes  vary  between  -0.54%  and  +1.7%.  In  two  cases, 
the  best  fit  straight  line  through  the  data  had  a  negative  slope, 
indicating  a  small  decrease  in  the  number  of  bubbles  counted  over 
time  since  irradiation.  These  values  indicate  that  the  bubble 
counts  for  all  the  devices  varied  very  little  with  time.  Although 
the  bubble  counts  for  the  devices  did  for  the  most  part  increase 
with  time,  it  can  be  concluded  from  this  test  that  the  bubbles, 
when  formed,  are  generally  large  enough  to  be  counted  by  the 
reader.  However,  some  of  the  bubbles  formed  initially  will  be  too 
small  to  be  counted.  A  very  small  number  of  bubbles  will  sometimes 
be  initially  too  small  and  will  be  below  the  minimum  sensitivity  of 
the  reader.  The  increase  or  decrease  in  number  of  bubbles  can  also 
be  partially  accounted  for  by  statistical  fluctuations  and  inherent 
reader  characteristics,  which  are  discussed  in  Section  5.5. 

The  second  group  of  ten  dosimeters  was  irradiated  once  using 
a  1  curie  Pu-Be  source.  For  this  test,  the  concept  of  the  "percent 
change  from  horizontal"  was  again  used  to  analyze  the  acquired 
data.  Figure  5-11  graphically  provides  the  percent  change  from 
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First  Run 


FIGURE  5-9.  Bar  graph  representing  the  percent  change  from  the 

horizontal  for  ten  devices. 


horizontal  for  the  ten  detectors  irradiated  by  the  Pu-Be  source. 
The  results  of  this  test  were  similar  to  those  obtained  using  the 
Cf^®^  source.  The  relative  percent  increases  in  the  counts  are  very 
small,  with  all  values  below  one  percent.  This  data  further 
indicates  that  the  bubble  growth  has  only  a  negligible  effect  upon 
the  reader's  ability  to  count  accurately  the  number  of  bubbles 
initially  present  after  irradiation.  Furthermore,  it  is  important 
to  note  that  the  higher  energy  of  the  Pu-Be  neutron  source  did  not 
noticeably  affect  the  bubble  growth  process. 

The  third  set  of  four  dosimeters  was  irradiated  using  the  Cf^®^ 
source  and  the  data  were  evaluated  over  a  very  long  period  of  time 
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FIGURE  5-10.  Bar  graph  representing  the  percent  change  from  the 
horizontal  for  ten  devices. 


following  the  irradiation.  Similar  to  the  previous  groups  of 
dosimeters,  the  data  were  analyzed  by  determining  the  percent 
change  from  the  horizontal  for  each  of  the  four  devices.  From  Fig. 
5-12  it  is  evident  that  the  devices  had  a  very  negligible  increase 
in  response  for  the  time  period  (1880  minutes)  considered.  The 
percent  change  for  device  G  is  the  largest  and  is  less  than  0.003 
percent.  From  this  data,  it  can  be  concluded  that,  over  long 
periods  of  time,  the  response  does  not  significantly  change.  When 
the  results  of  this  test  are  compared  with  the  results  of  the  first 
and  second  group  of  detectors  irradiated,  some  interesting 
conclusions  can  be  drawn. 
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FIGURE  5-11.  Bar  graph  representing  the  percent  change  from  the 
horizontal  for  ten  devices. 


The  change  in  the  bubble  counts  during  the  first  twelve 
minutes  after  irradiation  were  much  more  dramatic  than  the  change 
in  counts  for  the  extended  run.  This  makes  sense  when  considering 
the  definition  of  "percent  change  from  the  horizontal"  and  the 
drastic  difference  in  time  intervals  over  which  the  changes  were 
observed.  It  can  be  concluded  that  although  the  bubble  counts  did 
slightly  increase  over  the  first  twelve  minutes  after  irradiation, 
the  rate  of  increase  did  not  remain  constant,  but  rather  decreased 
with  time.  This  observation  remains  consistent  with  the  bubble 
growth  histories  provided  in  Fig.  5-7.  The  change  in  indicated 
number  of  bubbles  over  the  much  longer  extended  run  is  fairly 
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FIGURE  5-12 


Bar  graph  representing  the  percent  changes  from 
horizontal  for  four  detectors. 


negligible  when  compared  to  initial  increase  in  response  during  the 
first  twelve  minutes  after  irradiation.  Furthermore,  the  increase 
in  response  due  to  bubble  growth  was  found  to  be  much  less  than 
that  due  to  the  statistical  fluctuations  discussed  in  Section  5.1. 
As  a  result  of  these  tests,  it  appears  that  the  bubble  growth 
characteristics  of  the  BD-IOOR  bubble  detectors  have  a  very  minimal 
effect  upon  the  BTI  BDR-Series  II  optical  reader's  ability  to 
accurately  count  the  response.  Additionally,  the  BD-IOOR  bubble 
detectors  can  be  read  confidently  soon  after  irradiation. 


5 . 3  LINEARITY  STUDY 


One  of  the  most  important  characteristics  of  an  accurate 
radiation  detector  is  that  its  measurements  have  a  linear 
relationship  with  the  received  dose  or  fluence.  For  example,  if  a 
device  with  a  sensitivity  of  50  bubbles/mrem  is  exposed  to  a  source 
with  a  dose  rate  of  1  mrem/minute,  the  response  counted  must  be 
consistent  with  the  dose  the  device  received.  For  a  two  minute 
exposure,  the  device  should  have  approximately  100  bubbles,  while 
for  a  four  minute  exposure,  the  device  should  register  twice  as 
many  bubbles.  A  device  linearity  test  verifies  whether  the 
detector  sensitivity  (bubbles/mrem)  remains  constant  with  a 
changing  dose.  The  next  test  investigated  the  linear 
characteristics  of  the  BD-IOOR  bubble  detector  by  exposing  it  to 
different  neutron  doses. 


EXPERIMENTAL  SET-UP 

The  linearity  study  was  conducted  utilizing  ten  different  BD- 
lOOR  detectors.  The  devices  were  each  irradiated  for  time  periods 
of  two,  three,  four,  and  six  minutes  using  the  Cf^®^  source 
previously  described.  These  irradiation  periods  were  chosen 
because  they  were  expected  to  produce  bubble  counts  ranging  from 
200  to  about  600  bubbles,  or  great  enough  to  provide  reliable 
statistics,  and  still  lying  within  the  advertised  dynamic  range  of 
the  optical  reader.  The  actual  experimental  setup  for  the 
linearity  study  was  identical  to  that  outlined  in  Fig.  5-1.  All 
irradiations  were  performed  at  standard  atmospheric  pressure  and 
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constant  room  temperature  of  22 ®C. 


DISCUSSION  OF  RESULTS 

In  this  study,  ten  different  BD-IOOR  dosimeters  were  each 
exposed  to  the  Cf^®^  source  for  four  different  irradiation  time 
periods.  Figure  5-13  provides  a  graphical  display  of  the  data 
acquired  for  device  G,  one  of  the  ten  devices.  In  this  figure,  the 
number  of  bubbles  for  each  of  the  four  runs  is  plotted  with  respect 
to  the  length  of  irradiation  in  minutes.  The  best  fit  line  through 
the  data  was  determined  by  performing  a  linear  regression. 

A  linear  regression  was  also  performed  on  the  data  for  each  of 
the  other  nine  devices.  The  linearity  of  each  device's  behavior 
was  ascertained  by  calculating  the  linear  correlation  factor.  This 
parameter  ranges  between  zero  and  one,  with  one  being  exactly 
linear  (all  data  points  on  the  same  straight  line) ,  and  provides  an 
indication  of  how  closely  the  accumulated  data  matches  the  curve¬ 
fitting  model  used.  Figure  5-14  shows  the  correlation  factors  for 
each  of  the  ten  devices  used  in  the  linearity  study.  From  the 
results  provided  in  Fig.  5-14,  it  can  be  concluded  that  the  BD-IOOR 
detectors  behave  linearly  with  varying  dose,  for  the  dose  range 
covered  by  these  tests. 


5.4  TEMPERATURE  RESPONSE  STUDY 

The  critical  amount  of  energy,  E^,  necessary  to  form  a  bubble 
within  a  device  depends  upon  the  temperature  of  the  superheated 
liquid  droplet.  This  phenomenon  was  described  in  Chapters  II,  III, 
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Device  G 


and  IV  and  is  clearly  demonstrated  in  Fig  3-1.  As  a  consequence, 
the  response  of  a  bubble  detector  is  usually  temperature  sensitive 
and  non-linear ly  increases  with  increasing  temperature.  Prior 
studies*' have  verified  that  both  the  Apfel  detectors  and  earlier 
models  of  the  BD-IOOR  do  adhere  to  the  temperature  dependence 
characteristics  of  bubble  detectors. 

BTI  has  conducted  extensive  research  into  ways  of  redesigning 
♦  the  BD-IOOR  so  that  its  temperature  dependence  is  either  eliminated 

or  limited  to  a  certain  temperature  region.  As  mentioned  in 
Chapter  III  of  this  report,  one  of  the  compensation  techniques  is 
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FIGURE  5’'14.  Correlation  factors  for  each  of  the  ten  devices 


evaluated  in  the  linearity  study.  A  correlation 
factor  of  1.0  indicates  perfect  linearity. 


mechanical,  as  used  in  the  BD-IOOR.  A  compensating  material  which 
expands  and  contracts  with  changing  temperatures  is  placed  within 
the  device  immediately  above  the  rigid  gel  matrix  which  contains 
the  superheated  liquid  droplets.  Increasing  temperatures  cause  the 
compensating  material  to  expand  and  exert  a  pressure  upon  the  gel 
matrix.  As  a  result,  the  increased  pressure  counter-balances  the 
effects  of  increased  temperature  upon  E^.  This  study  evaluated  the 
temperature  characteristics  of  the  new  compensated  BD-IOOR  by 
measuring  the  response  of  the  device  at  many  different 
temperatures . 
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EXPERIMENTAL  SET-UP 

Ten  different  BD-IOOR  detectors  were  each  irradiated  for  five 
minutes  at  nine  different  temperatures  ranging  from  9.5®C  to 
46.1®C.  The  irradiations  were  performed  with  the  devices  held  at 
the  different  temperatures  by  placing  them  in  a  Fisher  Scientific 
low  temperature  incubator.  The  Cf^“  source  was  positioned 
immediately  above  the  center  of  the  incubator  so  that  each  device 
would  receive  an  equal  dose.  At  least  three  hours  of  wait  time 
were  allowed  between  each  run  so  that  both  the  incubator  and  BD- 
lOOR  detectors  could  reach  thermal  equilibrium.  Three  hours  were 
sufficient  to  ensure  thermal  equilibrium  of  the  superheated  liquid 
droplet  material  in  the  dosimeters,  based  on  experimental  work  by 
Harper,  who  monitored  the  heatup  and  cooldown  dynamic  temperature 
response  of  these  detectors  using  thermocouples. 


DISCUSSION  OF  RESULTS 

The  results  of  the  BD-lOOR's  temperature  compensation 
evaluation  are  presented  in  Fig.  5-15.  In  this  figure,  the 
relative  responses  of  both  the  uncompensated  BD-IOOR  and  the 
temperature-compensated  BD-IOOR  are  presented.  The  manufacturer 
provided  the  data  for  the  uncompensated  BD— lOOR.  Both  sets  of 
responses  were  normalized  to  the  arbitrary  value  of  1.0  at  a 
temperature  of  1.0  at  20°C.  The  relative  response  provided  in  Fig. 
5-15  accounts  for  the  responses  of  all  ten  devices  and  provides  a 
good  indication  of  the  overall  temperature  characteristics  of  the 


FIGURE  5-15.  The  experimentally  determined  relative  response  of 

the  temperature  compensated  BD-IOOR  compared  to 
that  of  the  uncompensated  BD-IOOR. 

compensated  BD-IOOR  device. 

From  Fig.  5-15,  it  appears  that  the  compensated  BD-IOOR  is 
superior  in  performance  to  the  earlier  uncompensated  version,  when 
considering  response  versus  temperature.  The  data  obtained  did  for 
the  most  part  verify  the  performance  specifications  provided  by  the 
manufacturer,  as  shown  in  Fig.  5-15 (A).  The  measured  flat  response 
was  only  a  few  degrees  centigrade  short  of  that  predicted  by  the 
manufacturer.  From  this  test,  it  can  be  concluded  that  the 
compensated  BD-IOOR  detector  appears  to  have  a  region  of  flat 
response  with  respect  to  temperature  for  a  relatively  large 


FIGURE  5-15 (A) .  Comparison  between  uncompensated  and 

compensated  BD-IOOR,  as  purported  by  the 
manufacturer ,  BTI . 


temperature  band  of  about  22 ®C  to  34 ®C. 

5 . 5  READER  EVALUATION 

Because  of  the  obvious  limitations  to  counting  bubbles  by  eye 
(inconsistency,  time-consuming,  tedious,  ,  and  of  marginal  accuracy 
as  the  number  of  bubbles  increases) ,  a  computer-driven,  optical 
reader  system,  the  BDR-Series  II  Bubble  Reader,  is  typically 
employed  in  applications  where  more  than  several  dozen  bubbles  are 
formed,  where  many  detectors  are  being  read,  or  where  permanent 
records  of  results  are  required.  All  of  the  bubble  counting 
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performed  for  the  purposes  of  this  report  involved  the  use  of  the 
BTI  BDR-Series  II  optical  reader  system.  Although  the  performance 
of  the  earlier  Series  I  reader  has  already  been  evaluated^*  as 
satisfactory,  little  has  been  published  about  the  characteristics 
of  this  new  reader. 

The  optical  reader  is  composed  of  two  cameras  which 
simultaneously  take  pictures  of  different  sections  of  the  bubble 
detector.  The  cameras  are  linked  to  a  computer  assembly  which 
analyzes  the  photos  by  examining  the  individual  pixels  which  make 
up  the  images.  Using  an  algorithm,  the  computer  is  able  to  count 
the  number  of  bubbles  by  differentiating  between  shadows  and  other 
bubbles.  The  computer  provides  different  threshold  settings  that 
can  be  manipulated  to  adjust  the  sensitivity  of  the  optical  reader. 

Initially,  there  were  several  unknowns  associated  with  the 
operation  of  the  optical  reader.  For  example,  it  was  not  clear 
which  threshold  setting  should  be  used  for  certain  quantities  of 
bubbles.  Furthermore,  the  dynamic  range  and  error  associated  with 
device's  readings  were  also  unknown.  Consequently,  these 
parameters  were  studied  in  order  to  gain  an  understanding  of  the 
reader's  operational  capabilities,  and  strengthen  the  conclusions 
drawn  from  the  data  collected  during  the  previously  described 
tests. 

EXPERIMENTAL  SET-UP 

A  group  of  tests  was  designed  to  determine  the  optimum 
threshold  setting  (explained  below) ,  the  reader's  dynamic  range. 
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and  any  error  associated  with  the  reader's  counts.  In  the  first 
series  of  tests,  three  previously  irradiated  BD”100R  detectors, 
each  containing  different  numbers  of  bubbles,  were  evaluated  at 
varying  reader  thresholds.  Two  of  the  devices  contained  a  small 
enough  number  of  bubbles  (namely  47  bubbles  and  87  bubbles) ,  that 
the  actual  number  of  bubbles  could  be  verified  by  eye.  This 
confirmation  was  done  by  visually  counting  the  bubbles  in  the 
devices  using  the  displayed  pictures  provided  by  the  reader's  video 
camera.  An  exact  correlation  was  also  seen  between  the  bubbles  as 
they  appeared  on  the  video  monitor  and  physically  visible  in  the 
detector  itself.  The  third  device  however,  contained  too  many 
bubbles,  approximately  140,  and  could  not  be  visually  verified. 
The  three  devices  were  each  read  at  separate  threshold  settings 
ranging  from  l  to  80.  Furthermore,  thirty  independent  readings 
were  taken  for  each  device  at  each  of  the  detector  thresholds 
evaluated. 

The  second  series  of  tests  utilized  ten  different  BD~100R 
dosimeters.  They  were  each  irradiated  for  the  time  periods  of  two, 
four,  six,  eight,  nine,  and  eleven  minutes  using  the  Cf^“  source 
pr©viously  described.  These  irradiation  periods  were  chosen  so 
that  the  full  dynamic  range  of  the  reader  (i.e.  maximum  number  of 
bubbles  that  can  accurately  be  counted)  could  be  determined  and 
evaluated.  The  experimental  set-up  was  identical  to  that  outlined 
in  the  linearity  study  and  in  Fig.  5-1.  All  irradiations  were 
performed  at  standard  atmospheric  pressure  and  constant  room 
temperature  of  22°C. 
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DISCUSSION  OF  RESULTS 


The  goals  of  the  first  series  of  tests  were  to  determine  the 
optimxim  threshold  setting  for  different  numbers  of  bubbles  and  to 
quantify  the  error  associated  with  the  reader's  counts.  When  the 
reader  system  creates  a  binary  image,  it  does  so  by  looking  at  each 
pixel  in  the  grey  scale  image  and  determining  if  its  grey  level 
value  is  greater  or  less  than  the  threshold  value.  If  it  is  less, 
the  corresponding  pixel  will  be  black  and  if  it  is  greater,  the 
binary  pixel  will  be  white.  Therefore,  this  threshold  value  is 
used  to  separate  bubbles  from  the  background  image.  If  the  bubble 
sample  change  dramatically,  it  may  be  necessary  to  adjust  this 
value  to  produce  the  optimal  binary  image. Threshold  value  may 
be  increased  or  decreased,  in  accordance  with  the  directions 
outlined  in  section  9  of  the  user's  manual,  using  the  +  or  -  keys 
of  the  numeric  keypad  until  the  operator  determines  that  all 
bubbles  are  being  outlined  with  solid  black  boundaries.  It  is 
possible  to  note  whether  the  reader  has  registered  a  particular 
bubble  or  not  by  observing  whether  a  small  white  square  appears  in 
the  center  of  the  bubble's  image  on  the  video  monitor.  Adjustments 
by  the  operator  are  intended  to  ensure  that  all  bona  fide  bubbles 
actually  get  counted,  while  non-bubbles  do  not  get  counted.  Figure 
5-16  shows  the  number  of  bubbles  counted  at  the  different  reader 
threshold  settings  for  the  device  containing  47  bubbles.  The 
number  of  bubbles  plotted  represents  the  average  of  the  thirty 
readings  counted  for  each  of  the  different  threshold  settings.  It 
is  important  to  note  that  between  the  threshold  settings  of  55  and 
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65,  a  plateau  of  relatively  constant  response  existed.  Furthermore, 
in  this  plateau,  the  reader  was  accurately  able  to  determine  the 
true  number  of  bubbles,  as  compared  to  counting  visually.  However, 
the  data  plotted  in  Fig.  5-16 
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FIGURE  5-16. 


The  average  number  of  bubbles  counted  for  various 
reader  threshold  settings.  The  detector  had 
exactly  47  bubbles. 


only  represents  the  average  of  thirty  independent  readings  taken  at 
each  threshold  setting.  The  amount  of  deviation  present  in  the 
data  was  also  analyzed  and  these  results  are  shown  in  Fig.  5—17. 
The  percent  standard  deviation  achieved  a  minimum  value  between  the 
threshold  settings  of  ,55  and  60.  For  these  same  threshold 
settings,  the  reader,  as  shown  in  Fig.  5-16,  was  found  to  be  the 


93 


most  accurate.  Thus  precision  and  accuracy  appear  to  be  optimal  at 


a  threshold  setting  between  55  and  60. 


FIGURE  5-17 .  The  percent  standard  deviation  of  the  counts  for 


various  reader  threshold  settings.  The  detector 
had  exactly  47  bubbles. 


Figure  5-18,  similar  to  Fig.  5-16,  shows  the  average  number  of 
bubbles  counted  at  each  of  the  different  threshold  settings  for  the 
device  containing  87  bubbles.  In  this  plot,  a  plateau  of 
relatively  constant  and  accurate  response  exists  between  the 
threshold  settings  of  42  and  62.  It  appears,  therefore,  that  the 
reader  was  able  to  accurately  detect  the  true  number  of  bubbles 
over  this  region  of  threshold  settings.  In  Fig.  5-19,  the  percent 
standard  deviation  of  the  counts  at  each  threshold  setting  is 
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presented.  The  unique  relationship  between  detector  accuracy  and 
the  minimum  standard  deviation  in  the  counts  was  observed  once 
again,  as  the  standard  deviation  of  the  counts  appeared  to  be 
minimal  at  a  threshold  setting  of  55  to  60,  which  is  within  the 
most  accurate  band  of  42  to  62. 
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FIGURE  5-18. 


The  average  number  of  bubbles  counted  for  various 
reader  threshold  settings.  The  detector  had  exactly 
87  bubbles. 


Figure  5-20  shows  the  average  counts  acquired  for  the  device 
containing  approximately  140  bubbles.  Although  the  actual  number 
of  bubbles  could  not  be  visually  verified  because  of  excessive 
bubble  density,  and  thus  compared  to  detector  performance,  a 
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plateau  of  consistent  response  did  develop  between  the  threshold 
settings  of  60  and  63.  This  plateau  did  correspond  to  a  bubble 
count  of  approximately  140  bubbles,  indicating  apparently  accurate 
reader  measurements.  Similar  to  the  previous  tests,  the  plot  of 
the  percent  standard  deviation  of  the  counts.  Fig.  5-21,  developed 
a  region  of  minimum  standard  deviation  (maximum  precision)  between 
the  settings  of  60  and  63,  which  corresponded  to  the  bubble  count 
plateau  shown  in  Fig.  5-20. 

As  a  result  of  the  first  series  of  tests,  several  conclusions 
can  be  drawn  about  the  operating  characteristics  of  the  optical 
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FIGURE  5-20.  The  average  number  of  bubbles  counted  for  various 

reader  threshold  settings .  The  detector  had 
approximately  140  bubbles. 

reader.  First  of  all,  different  threshold  settings  are  required 
for  different  numbers  of  bubbles.  This  is  not  inconsistent  with 
the  reader  system's  operating  instructions,  which  recommend 
adjusting  the  threshold  settings  for  bubbles  having  different 
•'bubble  characteristics"  or  for  detectors  with  vastly  different 
numbers  of  bubbles.  The  results  of  the  three  tests  indicate  that 
separate  "best"  threshold  settings  do  exist,  depending  on  the 
bubble  density,  for  which  the  reader's  accuracy  and  precision  is 
the  highest.  Next,  the  region  of  maximum  accuracy  corresponds 
directly  to  the  zone  in  which  the  smallest  percent  standard 
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FIGURE  5'*21.  The  percent  standard  deviation  of  the  counts  for 
various  reader  threshold  settings.  The  detector 
had  approximately  140  bubbles. 


deviation  in  the  counts  exists.  Essentially,  the  detector  becomes 
more  accurate  as  the  deviation  in  the  counts  decreases . 

The  second  series  of  tests  attempted  to  quantify  the  dynamic 
operating  range  of  the  reader  by  determining  its  upper  bounds. 
This  was  accomplished  by  repeating  the  linearity  study  and 
including  longer  periods  of  irradiation.  The  dynamic  limit  of  the 
reader  could  easily  be  determined  by  evaluating  the  number  of 
bubbles  at  which  the  detector  would  no  longer  provide  linear 
readings.  Figure  5-22  provides  the  results  of  the  second  linearity 
test.  In  this  plot,  the  average  counts  of  the  ten  different 
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FIGURE  5*22.  The  dynamic  limit  of  the  optical  reader  was  shown 
to  be  675  bubbles  for  ten  different  detectors. 


detectors  are  plotted  with  respect  to  the  length  of  irradiation. 
After  a  period  of  irradiation  of  nine  minutes,  the  devices  no 
longer  behaved  linearly.  For  this  test,  an  irradiation  time  of 
nine  minutes  corresponded  to  an  average  bubble  count  of  675 
bubbles.  Before  the  devices  contained  about  675  bubbles,  the 
reader  provided  consistently  linear  results,  as  demonstrated  by  a 
correlation  factor  of  0.9990.  Figure  5-23  emphasizes  the  linear 
behavior  of  the  devices  in  the  dynamic  region  by  showing  the 
correlation  factor  for  each  device. 
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Correlation  Factor 


FIGURE  5—23.  The  correlation  factors  for  each  of  the  ten  devices 


while  operating  in  the  dynamic  region  of  the 
optical  reader. 
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CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 

Some  overall  project  conclusions  are  presented  and  discussed. 
Areas  that  require  future  work  are  outlined. 

6 . 1  CONCLUSIONS 

At  the  outset  of  this  project,  several  problem  areas  were 
identified  relating  to  the  practical  operation  of  neutron  bubble 
detectors.  These  areas  included  the  strong  temperature  dependence 
associated  with  the  responses  of  the  Apfel  and  BD~100R  devices,  the 
concept  of  bubble  growth  and  how  it  relates  to  accurate  bubble 
counting,  and  also  the  low  sensitivities  inherently  associated  with 
the  devices.  The  report  was  begun  as  a  feasibility  study  which 
would  study  the  full  scope  and  possible  solutions  to  these 
problems . 

The  temperature  dependence  problem  was  viewed  as  perhaps  the 
most  damaging  disadvantage  associated  with  the  bubble  detectors. 
Consequently,  considerable  theoretical  work  was  undertaken  to  solve 
this  problem  and  redesign  the  devices  so  that  they  would  provide 
regions  of  consistent  response  with  changing  temperature.  As  a 
result  of  this  work,  a  protocol  was  established  for  determining 
specifically  whether  a  superheated  liquid  would  provide  a  flat 
response  over  a  specific  temperature  range.  Furthermore,  several 
alternate  compounds  were  identified  which  would  theoretically 
provide  a  flat  response.  Although  the  temperature  compensated  BD- 
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lOOR  detector  was  experimentally  shown  to  provide  a  plateau  of 
consistent  response,  the  concept  of  changing  the  droplet  material 
is  more  favorable.  Changing  the  droplet  material  solves  the  root 
of  the  problem,  and  the  compensation  does  not  rely  upon  mechanical 
methods  which  are  subject  to  failure.  Another  advantage  to  varying 
the  droplet  materials  is  that  by  doing  so,  the  devices  can  be 
designed  to  meet  the  specific  temperature  requirements  of  the  user, 
within  limits. 

A  series  of  experiments  were  conducted  at  the  United  States 
Naval  Academy's  neutron  generation  facility  to  test  the  suitability 
of  the  devices  to  measure  neutrons.  The  experiments  conducted 
include  the  repeatability  study,  the  bubble  growth  study,  the 
linearity  study,  the  study  of  the  new  temperature  compensated  BD- 
10 OR  detector,  and  the  reader  evaluation  study.  The  results  of 
these  tests  were  very  favorable,  indicating  that  the  BD-IOOR 
dosimeter  is  a  very  capable  neutron  detection  device.  Linearity 
and  repeatability  were  both  demonstrated,  within  the  limits  of  the 
experimental  parameters.  Furthermore,  the  bubble  growth  study 
demonstrated  that  although  the  bubbles  do  continue  to  grow  after 
their  initial  formation,  the  growth  has  a  marginal  effect  upon  the 
optical  reader's  ability  to  measure  accurately  the  number  of 
bubbles  in  the  detector  tubes.  An  in-depth  evaluation  of  the 
optical  reader  system  indicated  that,  although  adequate,  the 
operating  instructions  given  in  the  user's  manual  are  nebulous  in 
some  areas,  including  the  use  of  the  threshold  setting.  There  was 
also  no  mechanism  to  compare  the  optical  reader  system  output 
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(number  of  bubbles)  with  a  prime  standard  detector  of  known,  fixed 
number  of  bubbles,  thus  providing  an  alignment  of  the  system. 
There  was  also  no  provision  for  checking  the  calibration  between 
prime  standard  alignments,  or  during  extended  periods  of  use. 

Although  there  are  several  other  areas  that  must  be 
investigated  further,  preliminary  studies  indicated  that  the  bubble 
detector,  specifically  the  BD-IOOR,  is  a  valid  neutron  detector  and 
has  some  potential  to  be  used  as  a  detector  for  arms  control 
applications. 

6 . 2  FUTURE  WORK 

As  stated  in  Section  6.1,  there  are  several  other  topics 
concerning  detector  operation  that  must  be  explored  before  a  final 
decision  can  be  reached  concerning  the  suitability  of  the  BD^lOOR 
for  use  by  on-site  inspection  teams. 

One  of  the  primary  tasks  of  any  future  projects  will  be  to 
have  commercial  vendors  manufacture,  if  feasible,  new  bubble 
detectors.  These  detectors  would  incorporate  the  materials 
identified  in  this  report  as  having  potentially  superior 
temperature  response  characteristics.  Furthermore,  these  devices 
should  be  fully  tested  to  validate  previous  predictions  made 
concerning  their  temperature  response  characteristics  and  to 
evaluate  their  ability  to  measure  neutron  radiation  effectively. 
Preliminary  discussions  with  the  commercial  vendors,  specifically 
Apfel  Enterprises  and  BTI,  have  indicated  that  the  new  devices  can 
be  made . 
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Perhaps  the  most  important  task,  however,  is  to  develop 
operational  procedures  for  making  field  measurements  to  determine 
the  presence  or  absence  of  nuclear  re-entry  vehicles  contained 
within  a  weapon.  Once  the  procedures  have  been  determined,  they 
should  be  evaluated  using  a  laboratory  mockup  of  a  warhead  capable 
of  simulating  a  variety  of  configurations. 

Some  other  areas  of  investigation  should  include  a  complete 
evaluation  of  the  BD-IOOR'S  gamma  sensitivity.  Earlier  models  of 
the  device  were  found  to  be  gamma  insensitive,  however  recent 
changes  in  the  manufacturing  could  have  altered  this  insensitivity. 
Finally,  devices  with  increased  neutron  sensitivity  should  be 
ordered  from  the  vendors  in  order  to  evaluate  their 
characteristics.  Device  sensitivity  is  an  important  topic  because 
it  directly  affects  the  number  of  bubbles  that  a  device  will 
register  for  a  given  neutron  fluence.  It  is  favorable  to  have  the 
number  of  bubbles  received  for  a  given  dose  be  as  high  as  possible 
so  as  to  maximize  the  statistical  certainty  associated  with  the 
counts,  and  also  to  reduce  the  time  duration  of  an  exposure. 
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APPENDIX  A 


4 


Critical  Bubble  Data  for  Alternate  Compounds 

Data  is  graphically  presented  for  some  of  the  other  candidates 
studied  as  suitable  alternatives  to  the  CFC  type  superheated  liquid 
droplet  materials  contained  in  the  present  varieties  of  bubble 
neutron  dosimeters.  Apfel's  8DDIOO  ciirrently  employs  Freon-12  as 
the  superheated  liquid  droplet  material/  and  the  compound  utilized 
in  BTI's  BD-IOOR  is  proprietary,  and  unknown  at  this  time. 
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FIGURE  A-1. 


Critical  energy  and  bubble  radius  for  Freon-12  as  a 
function  of  temperature  at  1  atmosphere. 
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FIGURE  A-2.  Critical  energy  and  bubble  radius  for  Freon-13  as  a 

function  of  temperature  at  1  atmosphere. 
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FIGURE  A-3.  Critical  energy  and  bubble  radius  for  Freon-22  as  a 

function  of  temperature  at  1  atmosphere. 
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FIGURE  A-4 


Critical  energy  and  bubble  radius  for  HCFC-124  as  a 
function  of  temperature  at  1  atmosphere. 
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FIGURE  A-5.  Critical  energy  and  bubble  radius  for  Freon-114  as 

a  function  of  temperature  at  1  atmosphere. 


FIGURE  A-6. 


Critical  energy  and  bubble  radius  for  butane  as  a 
function  of  temperature  at  1  atmosphere. 


FIGURE  A-7.  Critical  energy  and  bubble  radius  for  HFC-125  as  a 

function  of  temperature  at  1  atmosphere. 
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FIGURE  A-8.  Critical  energy  and  bubble  radius  for  Freon-C318  as 
a  function  of  temperature  at  1  atmosphere. 
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